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Abstract 

In these lecture notes for the Les Houches School on Applications of Random Matri- 
ces in Physics we give an introduction to the connections between matrix models and 
topological strings. We first review some basic results of matrix model technology and 
then we focus on type B topological strings. We present the main results of Dijkgraaf 
and Vafa describing the spacetime string dynamics on certain Calabi-Yau backgrounds 
in terms of matrix models, and we emphasize the connection to geometric transitions 
and to large N gauge/string duality. We also use matrix model technology to analyze 
large N Chern-Simons theory and the Gopakumar-Vafa transition. 
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1 Introduction 

Topological string theory was introduced by Witten in H2j as a simplified model 
of string theory which captures topological information of the target space, and it has 
been intensively studied since then. There are three important lessons that have been 
learned in the last few years about topological strings: 

1) Topological string amplitudes are deeply related to physical amplitudes of type II 
string theory. 

2) The spacetime description of open topological strings in terms of string field theory 
reduces in some cases to very simple gauge theories. 

3) There is an open/closed topological string duality which relates open and closed 
string backgrounds in a precise way 

In these lectures we will analyze a particular class of topological string theories where 
the gauge theory description in (2) above reduces in fact to a matrix model. This was 
found by Dijkgraaf and Vafa in a series of beautiful papers [23 EHl 121], where they 
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also showed that, thanks to the connection to physical strings mentioned in (1), the 
computation of nonperturbative superpotentials in a wide class of Af = 1 gauge theories 
reduces to perturbative computations in a matrix model. This aspect of the work of 
Dijkgraaf and Vafa was very much explored and exploited, and rederived in the context 
of supersymmetric gauge theories without using the connection to topological strings. 
In these lectures we will focus on the contrary on (2) and (3), emphasizing the string 
field theory construction and the open/closed string duality. The applications of the 
results of Dijkgraaf and Vafa to supersymmetric gauge theories have been developed in 
many papers and reviewed for example in [Bj, and we will not cover them here. Before 
presenting the relation between matrix models and topological strings, it is worthwhile 
to give a detailed conceptual discussion of the general ideas behind (2) and (3) and 
their connections to large N dualities. 

In closed string theory we study maps from a Riemann surface S 9 to a target manifold 
X, and the quantities we want to compute are the free energies at genus g, denoted by 
F g (ti). Here, the U are geometric data of the target space X, and the free energies are 
computed as correlation functions of a two-dimensional conformal field theory coupled 
to gravity. In topological string theory there are two different models, the A and the 
B model, the target space is a Calabi-Yau manifold (although this condition can be 
relaxed in the A model), and the parameters tj are Kahler and complex parameters, 
respectively. The free energies are assembled together into a generating functional 

oo 

F(g s ,U) = Y.9 2 s 9 ~ 2F M> ( L1 ) 

9=0 

where g s is the string coupling constant. 

In open string theory we study maps from an open Riemann surface E 9i h to a target 
X, and we have to provide boundary conditions as well. For example, we can impose 
Dirichlet conditions by using a submanifold 5* of X where the open strings have to end. 
In addition, we can use Chan-Paton factors to introduce a U(N) gauge symmetry. The 
open string amplitudes are now F g> h, and in the cases that will be studied in these 
lectures the generating functional will have the form 

oo oo 

F(g s , N) = J2Y1 F g , h g 2 r 2+h N h . (1-2) 

9=0 h=l 

Physically, the introduction of Chan-Paton factors and boundary conditions through a 
submanifold S of X means that we are wrapping N (topological) D-branes around S. A 
slightly more general situation arises when there are n submanifolds Si, ■ ■ ■ ,S n where 
the strings can end. In this case, the open string amplitude is of the form F 9i h u ... } h n 
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and the total free energy is now given by 



oo 



oo 



F(g s ,N l )=J2 E F i 



g,hi>— ,h n 9 s 



(1.3) 



9=0 hi,— ,h n =l 



where h = Y17=i m ^ ne case °f °P en strings one can in some situations use string 
field theory to describe the spacetime dynamics. The open string field theory of Witten 
which was originally constructed for the open bosonic string theory, can also be 
applied to topological string theory, and on some particular Calabi-Yau backgrounds 
the full string field theory of the topological string reduces to a simple U(N) gauge 
theory, where g s plays the role of the gauge coupling constant and iV is the rank of 
the gauge group. In particular, the string field reduces in this case to a finite number 
of gauge fields. As a consequence of this, the open string theory amplitude F g h can 
be computed from the gauge theory by doing perturbation theory in the double line 
notation of 't Hooft [HE]- More precisely, F 9jh is the contribution of the fatgraphs of 
genus g and h holes. The idea that fatgraphs of a U(N) gauge theory correspond 
to open string amplitudes is an old one, and it is very satisfying to find a concrete 
realization of this idea in the context of a string field theory description of topological 
strings, albeit for rather simple gauge theories. 

The surprising fact that the full string field theory is described by a simple gauge 
theory is typical of topological string theory, and does not hold for conventional string 
models. There are two examples where this description has been worked out: 

1) The A model on a Calabi-Yau of the form X = T*M, where M is a three-manifold, 
and there are N topological D-branes wrapping M. In this case, the gauge theory is 
Chern-Simons theory on M |74j . 

2) The B model on a Calabi-Yau manifold X which is the small resolution of a 
singularity characterized by the hyperelliptic curve y 2 = (W'(x)) 2 . If W'(x) has degree 
n, the small resolution produces n two-spheres, and one can wrap Ni topological D- 
branes around each two-sphere, with i — 1, • • • , n. In this case Dijkgraaf and Vafa 
showed that the gauge theory is a multicut matrix model with potential W(x) \27\ . 

In both examples, the open string amplitudes F 9)h are just numbers computed by 
the fatgraphs of the corresponding gauge theories. 

The fatgraph expansion of a U(N) gauge theory can be resummed formally by in- 
troducing the so called 't Hooft parameter t = g s N. For example, in the case of the 
free energy, we can rewrite (jl.2J) in the form by defining 



oo 




(1.4) 



h=l 
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In other words, starting from an open string theory expansion we can obtain a closed 
string theory expansion by resumming the hole expansion as indicated in (|1.4j) . This 
idea was proposed by 't Hooft [66J and gives a closed string theory interpretation of a 
gauge theory. 

What is the interpretation of the above resummation for the gauge theories that de- 
scribe the spacetime dynamics of topological open string theories? As was explained in 
[35] (for the A model example above) and in ^3] (for the B model example), there is a 
geometric or large N transition that relates the open string Calabi-Yau background X 
underlying the gauge theory to a closed string Calabi-Yau background X' . The geomet- 
ric transition typically relates two different ways of smoothing out a singular geometry 
(the "resolved" geometry and the "deformed" geometry). Moreover, the "master field" 
that describes the large N limit jHH] turns out to encode the target space geometry of 
the closed string background, and the 't Hooft parameter becomes a geometric param- 
eter of the resulting closed geometry. The idea that an open string background with 
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Figure 1: This diagram summarizes the different relations between closed topological 
strings, open topological strings, and gauge theories. 

D-branes is equivalent to a different, geometric closed string background (therefore 
with no D-branes) appeared originally in the AdS/CFT correspondence In this 
correspondence, type IIB theory in flat space in the presence of D-branes is conjec- 
tured to be equivalent to type IIB theory in AdSs x S 5 with no D-branes, and where 
the radius of the S 5 is related to the 't Hooft parameter. The reason this holds is 
that, at large N, the presence of the D-branes can be traded for a deformation of the 
background geometry. In other words, we can make the branes disappear if we change 



5 



the background geometry at the same time. Therefore, as emphasized by Gopakumar 
and Vafa in (33] , large N dualities relating open and closed strings should be associated 
to transitions in the geometry. The logical structure of all the connections we have 
sketched is depicted in Fig. [T] 

In these lectures we will mostly focus on the B topological string, the Dijkgraaf-Vafa 
scenario, and the geometric transition of ^H]- For a detailed review of a similar story 
for the A string, we refer the reader to [HI]. The organization of these lectures is as 
follows. In section 2 we review some basic ingredients of matrix models, including 
saddle-point techniques and orthogonal polynomials. In section 3 we explain in detail 
the connection between matrix models and topological strings due to Dijkgraaf and 
Vafa. We first review the topological B model and its string field theory description, 
and we show that in the Calabi-Yau background associated to the resolution of a 
polynomial singularity, the string field theory reduces to a matrix model. We develop 
some further matrix model technology to understand all the details of this description, 
and we make the connection with geometric transitions. In section 4 we briefly consider 
the geometric transition of Gopakumar and Vafa from the point of view of the 
matrix model description of Chern-Simons theory. This allows us to use matrix model 
technology to derive some of the results of j3~H] . 

2 Matrix models 

In this section we develop some aspects and techniques of matrix models which will 
be needed in the following. There are excellent reviews of this material, such as for 
example pTll2^] . 

2.1 Basics of matrix models 

Matrix models are the simplest examples of quantum gauge theories, namely, they are 
quantum gauge theories in zero dimensions. The basic field is a Hermitian N x N 
matrix M. We will consider an action for M of the form: 



where g s and g p are coupling constants. This action has the obvious gauge symmetry 




(2.1) 



M -> UMU\ 

where U is a U(N) matrix. The partition function of the theory is given by 



(2.2) 




(2.3) 
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where the factor vo\(U(N)) is the usual volume factor of the gauge group that arises 
after fixing the gauge. In other words, we are considering here a gauged matrix model. 
The measure in the "path integral" is the Haar measure 



N 

N(N-i) -i— r i— r 

dM = 2-^—[[dM ii [I dReMijdlmMij. (2.4) 

i=l l<i<j<N 



The numerical factor in (J2.4)) is introduced to obtain a convenient normalization. 

A particularly simple example is the Gaussian matrix model, defined by the partition 
function 

Zg = 1 f d Me-W^ M \ (2.5) 



vo\(U(N)) 



We will denote by 



/ dMf(M) 



e 



-TrM 2 /2 9s 



the normalized vevs of a gauge-invariant functional f{M) in the Gaussian matrix 
model. This model is of course exactly solvable, and the vevs (J2.6|) can be computed 
systematically as follows. Any gauge-invariant function f(M) can be written as a linear 
combination of traces of M in arbitrary representations R of U(N). If we represent R 
by a Young tableau with rows of lengths Aj, with Ai > A 2 > • • ■ , and with £(R) boxes 
in total, we define the set of £(R) integers fi as follows 

fi = Xi + e{R)-i, i = l,-.. ,£(R). (2.7) 

Following [23], we will say that the Young tableau associated to R is even if the number 
of odd fi's is the same as the number of even /j's. Otherwise, we will say that it is 
odd. If R is even, one has the following result 



(Tr R M) G = c{R) dim R, (2.8) 

where 

c(R) = (-1)^ ^!^ (2.9) 

11/ odd,/' even 1/ / ) 

and A = £(R)/2 (notice that £{R) has to be even in order to have a non- vanishing 
result). Here dimi? is the dimension of the irreducible representation of SU(N) as- 
sociated to R, and can be computed for example by using the hook formula. On the 
other hand, if R is odd, the above vev vanishes. 

The partition function Z of more general matrix models with action (|2.1jl can be 
evaluated by doing perturbation theory around the Gaussian point: one expands the 
exponential of J2 P >3(.9p/ ' g s )^M p j p in (|2.Hjl . and computes the partition function as a 
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power series in the coupling constants g p . The evaluation of each term of the series 
involves the computation of vevs like ([2.60 . Of course, this computation can be inter- 
preted in terms of Feynman diagrams, and as usual the perturbative expansion of the 
free energy 



will only involve connected vacuum bubbles. 

Since we are dealing with a quantum theory of a field in the adjoint representation 
we can reexpress the perturbative expansion of F in terms of fatgraphs, by using the 
double line notation due to 't Hooft jHE]- The purpose of the fatgraph expansion is 
the following: in U(N) gauge theories there is, in addition to the coupling constants 
appearing in the model (like for example g s ,g p in (|2.ip ). a hidden variable, namely 
N, the rank of the gauge group. The N dependence in the perturbative expansion 
comes from the group factors associated to Feynman diagrams, but in general a single 
Feynman diagram gives rise to a polynomial in N involving different powers of N. 
Therefore, the standard Feynman diagrams, which are good in order to keep track of 
powers of the coupling constants, are not good in order to keep track of powers of N. 
If we want to keep track of the N dependence we have to "split" each diagram into 
different pieces which correspond to a definite power of N. To do that, one writes the 
Feynman diagrams of the theory as "fatgraphs" or double line graphs, as first indicated 
by 't Hooft [HE]- Let us explain this in some detail, taking the example of the matrix 



Figure 2: The index structure of the field in the adjoint representation of U(N) is 
represented through a double line. 

model with a cubic potential (i.e. g p = in (j2.1j) for p > 3). The fundamental field 
My is in the adjoint representation. Since the adjoint representation of U(N) is the 
tensor product of the fundamental N and the antifundamental N, we can look at % 
(resp. j) as an index of the fundamental (resp. antifundamental) representation. We 
will represent this double-index structure by a double line notation as shown in Fig. El 
The only thing we have to do now is to rewrite the Feynman rules of the theory by 
taking into account this double-line notation. For example, the kinetic term of the 
theory is of the form 



F = log Z 




(2.10) 
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This means that the propagator of the theory is 

(MyMu) = g s 8 u 5 jk (2.11) 
and can be represented in the double line notation as in Fig. |3J Next, we consider the 



k 6 ii 6 J k 



Figure 3: The propagator in the double line notation, 
vertices of the theory. For example, the trivalent vertex given by 



— Tr M 3 = — V Ma M jk M; 
9s 9s (jt 



hi 



(2.12) 



can be represented in the double line notation as in Fig. 0] A vertex of order p can 
be represented in a similar way by drawing p double lines joined together. Once 




Figure 4: The cubic vertex in the double line notation. 

we have rewritten the Feynman rules in the double-line notation, we can construct 
the corresponding graphs, which look like ribbons and are called ribbon graphs or 
fatgraphs. It is clear that in general a usual Feynman diagram can give rise to many 

different fatgraphs. Consider for example the one-loop vacuum diagram (3), which 
comes from contracting two cubic vertices. In the double line notation the contraction 
can be done in two different ways. The first one is illustrated in Fig. El and gives a 
factor 

J2 {M ij M mn ){M jk M pm ){M kl M np ) = g 3 s N 3 . (2.13) 

ijkmnp 
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Figure 5: Contracting two cubic vertices in the double line notation: the N 3 contribu- 
tion. 

The second one is shown in Fig. El and gives a factor 

{M ij M mn ){M jk M np ){M ki M pm ) = g 3 s N. (2.14) 

ijkmnp 

In this way we have split the original diagram into two different fat graphs with a well- 
defined power of N associated to them. The number of factors of N is simply equal to 
the number of closed loops in the graph: there are three closed lines in the fatgraph 
resulting from the contractions in Fig. |5] (see the first graph in Fig. 0), while there is 
only one in the diagram resulting from Fig. El In general, fatgraphs turn out to be 




Figure 6: Contracting two cubic vertices in the double line notation: the N contribu- 
tion. 

characterized topologically by the number of propagators or edges E, the number of 
vertices with p legs V p , and the number of closed loops h. The total number of vertices 
is V — Yip Vp- Each propagator gives a power of g s , while each interaction vertex with 
p legs gives a power of g p /g s - The fatgraph will then give a factor 

p 
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The key point now is to regard the fatgraph as a Riemann surface with holes, in which 
each closed loop represents the boundary of a hole. The genus g of such a surface is 
determined by the elementary topological relation 

2g-2 = E-V -h (2.16) 

therefore we can write (|2.15jl as 

g 2 r 2+h N h J] # = 9 2 r 2 t h II # (2-17) 
v v 

where we have introduced the 't Hooft parameter 

t = Ng s (2.18) 

The fatgraphs with g — are called planar, while the ones with g > are called 
nonplanar. The graph giving the N 3 contribution in Fig. |3] is planar: it has E = 3, 
V3 = 2 and ft = 3, therefore g = 0, and it is a sphere with three holes. The graph in 
Fig. EI is nonplanar: it has E — 3, V3 = 2 and ft, = 1, therefore g — 1, and represents a 
torus with one hole (it is easy to see this by drawing the diagram on the surface of a 
torus) . 

We can now organize the computation of the different quantities in the matrix model 
in terms of fatgraphs. For example, the computation of the free energy is given in the 
usual perturbative expansion by connected vacuum bubbles. When the vacuum bubbles 
are written in the double line notation, we find that the perturbative expansion of the 
free energy is given by 

00 00 

f = Y.H F ^ 2 r 2 t\ (2.19) 

g=0 h=l 

where the coefficients F gi h (which depend on the coupling constants of the model g p ) 
takes into account the symmetry factors of the different fatgraphs. We can now formally 
define the free energy at genus g, F g (t), by keeping g fixed and summing over all closed 
loops ft as in (jl.4[) . so that the total free energy can be written as 

00 

F = J2 F ^9s 9 ~ 2 - (2-20) 

9=0 

This is the genus expansion of the free energy of the matrix model. In ()2.20|) we have 
written the diagrammatic series as an expansion in g s around g s = 0, keeping the 't 
Hooft parameter t = g s N fixed. Equivalently, we can regard it as an expansion in 1/N, 
keeping t fixed, and then the N dependence appears as N 2 ~ 2g . Therefore, for t fixed 
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and N large, the leading contribution comes from planar diagrams with g = 0, which 
go like 0(N 2 ). The nonplanar diagrams give subleading corrections. Notice that F g (t), 
which is the contribution to F to a given order in g s , is given by an infinite series where 
we sum over all possible numbers of holes h, weighted by t h . 



Figure 7: Two planar diagrams in the cubic matrix model. 

Example. One can show that 

((Tr M 3 ) 2 ) G = g 3 s (12N 3 + 3N), 

where the first term corresponds to the two planar diagrams shown in Fig. [7| (contribut- 
ing 3iV 3 and 9N 3 , respectively), and the second term corresponds to the nonplanar 
diagram shown in Fig. |H1 Therefore, in the cubic matrix model the expansion of the 
free energy reads, at leading order, 

F-F G = ^g s g 2 N 3 + +■■■ (2.21) 

There is an alternative way of writing the matrix model partition function which is 
very useful. The original matrix model variable has N 2 real parameters, but using the 
gauge symmetry we can see that, after modding out by gauge transformations, there 
are only N parameters left. We can for example take advantage of our gauge freedom 
to diagonalize the matrix M 

M -> UMU ] = D, (2.22) 

with D = diag(Ai, • • ■ , Ajv), impose this as a gauge choice, and use standard Faddeev- 
Popov techniques in order to compute the gauge- fixed integral (see for example 0). 
The gauge fixing (|2.22|) leads to the delta-function constraint 

8( v M)=Y[8®( u M ij ) (2.23) 

i<3 
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where U M = UMW. We then introduce 

A~ 2 (M) = / dU5( u M). (2.24) 



It then follows that the integral of any gauge-invariant function f(M) can be written 

as 

N 

J dMf(M) = J dMf(M)A 2 (M) J dUS( u M) = tt N J f[ dX t A 2 (X)f(X), (2.25) 

where we have used the gauge invariance of A(M), and 

Q N = dU (2.26) 



is proportional to the volume of the gauge group U(N), as we will see shortly. We have 
to evaluate the the factor A(A), which can be obtained from (|2.24|) by choosing M to 
be diagonal. If 

F(M) = 

is the gauge-fixing condition, the standard Faddeev-Popov formula gives 



/ 5F( U M) 



A 2 (M)=det(^A_- j (2.27) 

F=0 



where we write U = e A , and A is a anti-Hermitian matrix. Since 

Fij( u D) = (UDU% = AijiXi — Xj) -\ • (2.28) 

()2.27|) leads immediately to 

A 2 (A) = H(X t - Xj) 2 , (2.29) 



i<j 



the square of the Vandermonde determinant. Finally, we fix the factor a s follows. 
The Gaussian matrix integral can be computed explicitly by using the Haar measure 
(j2.4)l . and is simply 

J dMe~^ M ' 2 = (27rg s ) N2 / 2 . (2.30) 
On the other hand, by ()2.25|) this should equal 

r N 

Vl N / Y[dXiA 2 (X)e-^^ x l (2.31) 

J i=l 
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The integral over eigenvalues can be evaluated in various ways, using for example the 
Selberg function (HSj or the technique of orthogonal polynomials that we describe in 
the next subsection, and its value is 

gf / 2 (2vr) 7V / 2 G 2 (iV + 2) (2.32) 

where G 2 (z) is the Barnes function, defined by 

G 2 (z + l)=T(z)G 2 (z), G 2 (l) = l. (2.33) 

Comparing these results, we find that 

N(N-1) 

n N = {2n 2 (2.34) 

Using now (see for example [fiTi]): 

"™ = WnW (2 ' 35) 

we see that 

J dMf(M) = ±^ / n<O.A'(A)/(A). (2.36) 

The factor N\ in the r.h.s. of (|2.36|) has an obvious interpretation: after fixing the gauge 
symmetry of the matrix integral by fixing the diagonal gauge, there is still a residual 
symmetry given by the Weyl symmetry of U(N), which is the symmetric group Sn 
acting as permutation of the eigenvalues. The "volume" of this discrete gauge group is 
just its order, \Sn\ = N\, and since we are considering gauged matrix models we have 
to divide by it as shown in (|2.36|) . As a particular case of the above formula, it follows 
that one can write the partition function ()2.3)1 as 

11 r N 

z = Mp^F / n^^w^'"'' 1 " < 2 - 37 > 

The partition function of the gauged Gaussian matrix model (|2.5J) is given essentially 
by the inverse of the volume factor. Its free energy to all orders can be computed by 
using the asymptotic expansion of the Barnes function 

TV" 2 1 3 1 

logG 2 (iV + l) = — log N-— log N--N 2 + -N log 2tt + C'(-l) 

00 

y & 2g 2 _ 2 
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where B 2g are the Bernoulli numbers. Therefore, we find the following expression for 
the total free energy: 



N 2 / 3\ 1 

1 -(log(iV^)-0--logiV + C / (-i; 



oo 

y V2 9 N 2- 2g _ ( } 



If we now put N = t/g s , we obtain the following expressions for F g (t): 

Fo(t) = l t2 { l °St-l), 
Fi(t) = ~logt, 

9K) 2g(2g-2) 

2.2 Matrix model technology I: saddle-point analysis 

The computation of the functions F g (t) in closed form seems a difficult task, since in 
perturbation theory they involve summing up an infinite number of fatgraphs (with 
different numbers of holes h). However, in the classic paper it was shown that, 
remarkably, Fo(t) can be obtained by solving a Riemann-Hilbert problem. In this 
section we will review this procedure. 

Let us consider a general matrix model with action W(M), and let us write the 
partition function after reduction to eigenvalues (|2.37|) as follows: 

„ 1 f^rdX 



J f[ ^_ e N*S cBW (2 40) 



i=l 

where the effective action is given by 

1 - 2 

s eff (A) = -— + jp E lo § i A * - A ii- ( 2 - 41 ) 

i=l i<j 

Notice that, since a sum over N eigenvalues is roughly of order N, the effective action 
is of order 0(1). We can now regard N 2 as a sort of fi~ x in such a way that, as iV — > oo, 
the integral ()2.40|) will be dominated by a saddle-point configuration that extremizes 
the effective action. Varying S e s(X) w.r.t. the eigenvalue Aj, we obtain the equation 
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The eigenvalue distribution is formally defined for finite N as 



1 - 

where the Aj solve (|2.42j) . In the large iV limit, it is reasonable to expect that this 
distribution becomes a continuous function with compact support. We will assume 
that p(X) vanishes outside an interval C. This is the so-called one-cut solution. 

Qualitatively, what is going on is the following. Assume for simplicity that W(x), 
the potential, has only one minimum . We can regard the eigenvalues as coordinates 
of a system of N classical particles moving on the real line. The equation ()2.42|) says 
that these particles are subject to an effective potential 

It 

W e z{\) = W(Xi) - - ^ log |Ai - \j\ (2.44) 

which involves a logarithmic Coulomb repulsion between eigenvalues. For small 't 
Hooft parameter, the potential term dominates over the Coulomb repulsion, and the 
particles tend to be in the minimum X* of the potential W'(x*) = 0. This means that, 
for t = 0, the interval C collapses to the point x*. As t grows, the Coulomb repulsion 
will force the eigenvalues to be apart from each other and to spread out over an interval 
C. 

We can now write the saddle-point equation in terms of continuum quantities, by 
using the standard rule 



/ fwp^ dX - ( 2 - 45 ) 

i=l ^ c 

Notice that the distribution of eigenvalues p(X) satisfies the normalization condition 

J p{X)dX = 1. (2.46) 

The equation ()2.42|) then becomes 

where P denotes the principal value of the integral. The above equation is an integral 
equation that allows one in principle to compute p(A), given the potential W(X), as a 
function of the 't Hooft parameter t and the coupling constants. Once p(A) is known, 
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one can easily compute F (t): in the saddle-point approximation, the free energy is 
given by 

±F = S cS (p) + 0(N- 2 ), (2.48) 
where the effective action in the continuum limit is a functional of p: 

S eS ( p ) = -- [ d\p{X)W{\)+ I dXdX'p(X)p(X')\og\X- A'|. (2.49) 

t JC JcxC 

Therefore, the planar free energy is given by 

F (t)=t 2 S cS (p), (2.50) 

Since the effective action is evaluated on the distribution of eigenvalues which solves 
()2.47|) . one can simplify the expression to 

Fo(t) = -\f d\p(X)W(X). (2.51) 

Similarly, averages in the matrix model can be computed in the planar limit as 

-jy{TrM e ) = j dXX l p{X). (2.52) 

We then see that the planar limit is characterized by a classical density of states p(X), 
and the planar piece of quantum averages can be computed as a moment of this density. 
The fact that the planar approximation to a quantum field theory can be regarded as a 
classical field configuration was pointed out in |nH| (see [2II| for a beautiful exposition). 
This classical configuration is often called the master field. In the case of matrix models, 
the master field configuration is given by the density of eigenvalues p(A), and as we 
will see later it can be encoded in a complex algebraic curve with a deep geometric 
meaning. 

The density of eigenvalues is obtained as a solution to the saddle-point equation 
(I2.47j) . This equation is a singular integral equation which has been studied in detail 
in other contexts of physics (see, for example, jSZ])- The way to solve it is to introduce 
an auxiliary function called the resolvent. The resolvent is defined as a correlator in 
the matrix model: 

w(p) = Tf (Tr j-=), (2.53) 

which is in fact a generating functional of the correlation functions ()2.52j) : 

k=0 
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Being a generating functional of connected correlators, it admits an expansion of the 
form |2D|: 

oo 

c^(p) = 5>> fl (p), (2-55) 

9=0 

and the genus zero piece can be written in terms of the eigenvalue density as 

MP) = I dX-^r (2.56) 
J P~ A 

The genus zero resolvent ()2.56|) has three important properties. First of all, as a 
function of p it is an analytic function on the whole complex plane except on the interval 
C, since if A G C one has a singularity at A = p. Second, due to the normalization 
property of the eigenvalue distribution (j2.46j) . it has the asymptotic behavior 

LU (p) ~ -, P -> oo- (2.57) 
P 

Finally, one can compute the discontinuity of u>o{p) as one crosses the interval C. This 
is just the residue at A = p, and one then finds the key equation 

p(A) = ( Wo (A + ie) - w„(A - ie)) . (2.58) 

27U 

Therefore, if the resolvent at genus zero is known, the eigenvalue distribution follows 
from (|2.58|) . and one can compute the planar free energy. On the other hand, by looking 
again at the resolvent as we approach the discontinuity, we see that the r.h.s. of ()2.47|) 
is given by — (uo(p + ie) + uo(p — ie))/2, and we then find the equation 

u (p + ie) + u Q (p - ie) = --W'(p), (2.59) 

which determines the resolvent in terms of the potential. In this way we have reduced 
the original problem of computing F (t) to the Riemann-Hilbert problem of computing 
o>o(A). There is in fact a closed expression for the resolvent in terms of a contour 
integral jSB] which is very useful. Let C be given by the interval b < A < a. Then, one 
has 

^q/iigfip-'y-;)) 1 (2 .60) 

w 2t J e 2m p-z \(z-a)(z-b) J v ; 

This equation is easily proved by converting ()2.59|) into a discontinuity equation: 

1 W'(p) 

uio(p + ie)-Q (p-ie) = -- , (2.61) 

t V(p-a)(p-6) 
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where ujq(p) = uj (p)/^(p — a) (p — b). This equation determines tu (p) to be given 
by (|2.6U|) up to regular terms, but because of the asymptotics (|2.57|) . these regular 
terms are absent. The asymptotics of Wq(p) also gives two more conditions. By taking 
p — > oo, one finds that the r.h.s. of (j2.6()j) behaves like c + d/p + 0{l/p 2 ). Requiring 
the asymptotic behavior (|2.57|) imposes c = and d = 1, and this leads to 



dz W'(z) 



2?ri y/( z -a)(z-b) 
dz zW'(z) 
2riy/{z-a){z-b) 



2t. (2.62) 



These equations are enough to determine the endpoints of the cuts, a and b, as functions 
of the 't Hooft coupling t and the coupling constants of the model. 

The above expressions are in fact valid for very general potentials (we will apply 
them to logarithmic potentials in section 4), but when W(z) is a polynomial, one can 
find a very convenient expression for the resolvent: if we deform the contour in ()2.60|) 
we pick up a pole at z — p, and another one at infinity, and we get 

"o(p) = ^W'(p) - ^(p-a)(p-b)M(p), (2.63) 

where 

„„, , f dz W'il/z) 1 , nnA . 

M(p) = <b ' 2.64 

J 2m 1-pz y/(l-az){l-bz) 

Here, the contour is around z = 0. These formulae, together with the expressions 
(12.62)1 for the endpoints of the cut, completely solve the one-matrix model with one 
cut in the planar limit, for polynomial potentials. 

Another way to find the resolvent is to start with (J2.42)) . multiply it by l/(Aj —p), 
and sum over i. One finds, in the limit of large N, 

(u (p)) 2 - jW'(p)u (p) + -^R(p) = 0, (2.65) 

where 

R(p) = At / d\p(X)—^ (2.66) 

J P- A 

Notice that ()2.65|) is a quadratic equation for uo(p) and has the solution 

w (p) = Yt { W '( p) ~ ^(^'(P)) 2 - R(P)) , (2-67) 
which is of course equivalent to ([2.63)1 . 
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A useful way to encode the solution to the matrix model is to define 



y(p) = W'(p)-2tuj (p). (2.68) 
Notice that the force on an eigenvalue is given by 

fip) = -KM = -\{V(P + ie) + vip - ie))- (2-69) 

In terms of y(p), the quadratic equation ()2.fi5j) determining the resolvent can be written 

as 

y 2 = W'(p) 2 -R(p). (2.70) 

This is nothing but the equation of a hyperelliptic curve given by a certain deformation 
(measured by Rip)) of the equation y 2 = W'ip) 2 typical of singularity theory. We 
will see in the next section that this result has a beautiful interpretation in terms of 
topological string theory on certain Calabi-Yau manifolds. 

Example. The Gaussian matrix model. Let us now apply this technology to the 
simplest case, the Gaussian model with W{M) = M 2 /2. Let us first look for the 
position of the endpoints from (|2.62j) . Deforming the contour to infinity and changing 
z — ► l/z, we find that the first equation in (|2.62j) becomes 

dz 1 1 . 

o ; = 0, 2.71 

2m z 2 y/(l-az){l-bz) 

where the contour is now around z = 0. Therefore a + b = 0, in accord with the 
symmetry of the potential. Taking this into account, the second equation becomes: 

dz 1 1 , _ . 
? , = 2t, 2.72 

2<ki z 3 VI - a 2 z 2 V ; 

and gives 

a = 2\Tt. (2.73) 

We see that the interval C = [—a, a] = [-2y/i, 2\Tt] opens as the 't Hooft parameter 
grows up, and as t — > it collapses to the minimum of the potential at the origin, as 
expected. We immediately find from (j2.6r>|) 

Mp) = Y t (p-v / P 2 - At )> ( 2 - 74 ) 

and from the discontinuity equation we derive the density of eigenvalues 

p(X) = ±-VAt^. (2.75) 
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The graph of this function is a semicircle of radius 2y/i, and the above eigenvalue 
distribution is the famous Wigner-Dyson semicircle law. Notice also that the equation 
(I2.70J1 is in this case 

y 2 =p 2 - At. (2.76) 

This is the equation for a curve of genus zero, which resolves the singularity y 2 = p 2 . 
We then see that the opening of the cut as we turn on the 't Hooft parameter can be 
interpreted as a deformation of a geometric singularity. This will be later interpreted in 
section 3.5 from the point of view of topological string theory on Calabi-Yau manifolds. 

Exercise. Resolvent for the cubic matrix model. Consider the cubic matrix model 
with potential W(M) = M 2 /2 + g 3 M 3 /3. Derive an expression for the endpoints of 
the one-cut solution as a function of t, g^, and find the resolvent and the planar free 
energy. The solution is worked out in [T2*j . 

Although we will not need it in this review, there are well-developed techniques to 
obtain the higher genus F g (t) as systematic corrections to the saddle-point result F (t) 
[SIIS21- Interestingly enough, these corrections can be computed in terms of integrals 
of differentials defined on the hyperelliptic curve ([2.70)1 . 

We have so far considered the so-called one cut solution to the one-matrix model. 
This is not, however, the most general solution, and we now will consider the multicut 
solution in the saddle-point approximation. Recall from our previous discussion that 
the cut appearing in the one-matrix model was centered around a minimum of the 
potential. If the potential has many minima, one can have a solution with various 
cuts, centered around the different minima. The most general solution has then s cuts 
(where s is lower or equal than the number of minima n), and the support of the 
eigenvalue distribution is a disjoint union of s intervals 

C = U s l=1 C h (2.77) 

where 

Ci = [x 2 i,x 2 i-i] (2.78) 

and %2 S < ■ ■ ■ < x\. The equation (|2.fi7j) still gives the solution for the resolvent, and it 
is easy to see that the way to have multiple cuts is to require Uo(p) to have 2s branch 
points corresponding to the roots of the polynomial W'(z) 2 — R(z). Therefore we have 



which can be solved in a compact way by 



2s 



* H(p-x k )M(p), (2.79) 
\ fc=i 



i r dzW(z) f-Xp-x k Y 
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In order to satisfy the asymptotics (|2.57|) the following conditions must hold: 

hs = -i— / W ' {Z) r , 1 = 0,1,-, s. (2.81) 

In contrast to the one-cut case, these are only s + 1 conditions for the 2s variables 
Xk representing the endpoints of the cut. For s > 1, there are not enough conditions 
to determine the solution of the model, and we need extra input to determine the 
positions of the endpoints Xk- Usually, the extra condition which is imposed is that 
the different cuts are at equipotential lines (see for example [HJI1]). It is easy to see 
that in general the effective potential is constant on each cut, 

W cS (p) = Ti, pe C h (2.82) 

but the values of will be in general different for the different cuts. This means that 
there can be eigenvalue tunneling from one cut to the other. The way to guarantee 
equilibrium is to choose the endpoints of the cuts in such a way that 1^ = F for all 
i = 1, ■ ■ ■ ,s. This gives the s — 1 conditions: 

W eS (x 2i+ i) = W eS (x 2i ), z = l,-..,s-l, (2.83) 

which, together with the s + 1 conditions ()2.81|) provide 2s constraints which allow one 
to find the positions of the 2s endpoints X{. We can also write the equation ()2.83j) as 

/ dzM(z) Y[(z-x k )^ = 0. (2.84) 

Jx3i+1 k= i 

In the context of the matrix models describing topological strings, the multicut solution 
is determined by a different set of conditions and will be described in section 3.4. 

2.3 Matrix model technology II: orthogonal polynomials 

Another useful technique to solve matrix models involves orthogonal polynomials. This 
technique was developed in [HI IHj (which we follow quite closely) , and provides explicit 
expressions for F g (t) at least for low genus. This technique turns out to be particularly 
useful in the study of the so-called double-scaling limit of matrix models |13J. We will 
use this technique to study Chern-Simons matrix models, in section 4, therefore this 
subsection can be skipped by the reader who is only interested in the conventional 
matrix models involved in the Dijkgraaf-Vafa approach. 

The starting point of the technique of orthogonal polynomials is the eigenvalue rep- 
resentation of the partition function 

' J 8=1 
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where W(X) is an arbitrary potential. If we regard 

df, = e-^ WW ^ (2.86) 
as a measure in R, one can introduce orthogonal polynomials p n (X) denned by 

dnp n (X)p m (X) = h n 5 nm} n > 0, (2.87) 



where p n {X) are normalized by requiring the behavior p n (X) = X n + ■ • ■ . One can now 
compute Z by noting that 

A(A) = detp i _i(A i ). (2.88) 
By expanding the determinant as 

^(-l) E W[I^H(A t ) (2.89) 

where the sum is over permutations o of iV indices and e(cr) is the signature of the 
permutation, we find 

N-l N 

Z=J[h i = hS]lr?- i , (2.90) 

i=0 i=l 

where we have introduced the coefficients 

r fc = ~r ~~ j fc > 1. (2.91) 

One of the most important properties of orthogonal polynomials is that they satisfy 
recursion relations of the form 

(A + s n )p n (X) = p n +i(X) + r nPn ^(X). (2.92) 

It is easy to see that the coefficients r n involved in this relation are indeed given by 
()2.91|) . This follows from the equality 

K+i = J dfJ,p n+ i(X)Xp n (X), (2.93) 

together with the use of the recursion relation for Xp n+ i(X). For even potentials, s n = 0. 

As an example of this technique, we can consider again the simple case of the Gaus- 
sian matrix model. The orthogonal polynomials of the Gaussian model are well-known: 
they are essentially the Hermite polynomials H n {x), which are defined by 

Jn 

H n (x) = (-l)-e^e-* 2 . (2.94) 
ax n 
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More precisely, one has 



Pn{X) 



n/2 



H n (x/J2g~ a ), 



and one can then check that 



h n = {^y^9s, r G n =ng s . 



(2.95) 



(2.96) 



Using now ()2.90|) we can confirm the result ()2.33|) that we stated before. 

It is clear that a detailed knowledge of the orthogonal polynomials allows the com- 
putation of the partition function of the matrix model. It is also easy to see that the 
computation of correlation functions also reduces to an evaluation in terms of the co- 
efficients in the recursion relation. To understand this point, it is useful to introduce 
the orthonormal polynomials 

1 



V n (X) 

which satisfy the recursion relation 



=Pn(A), 



\V n {\) = -S n V n (X) + 0VH7Pn+l(A) + V^ n -l(A). 



(2.97) 



(2.98) 



Let us now consider the normalized vev (Tr M e ) , which in terms of eigenvalues is given 
by the integral 



1 r N J\ N 

^^) = wb/n^ w(A ' ) ^ 2 w(E^- 

J 1=1 ' 1=1 



(2.99) 



By using (|2.88|) it is easy to see that this equals 



N-l 

£ / ^l(A). 

j=0 J 



;2.ioo) 



This integral can be computed in terms of the coefficients in (J2.97j) . For example, for 
t = 2 we find 

N-l 

(TTM 2 ) = Y,^ + r J+1 + r] ), (2.101) 

j=0 

where we put r = 0. A convenient way to encode this result is by introducing the 
Jacobi matrix 



J 



( 


1/2 





■ 


• A 







1/2 

r 2 


■ 







1/2 





1/2 




v ••■ 








■7 



;2.102) 
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as well as the diagonal matrix 



It then follows that 



S 



/s 
si ••• 
s 2 ••■ 



(2.103) 



(Tr M l ) = Tr (J - S) e . (2.104) 



The results we have presented so far give the exact answer for the correlators and the 
partition function, at all orders in 1/N. As we have seen, we are particularly interested 
in computing the functions F g (t) which are obtained by resumming the perturbative 
expansion at fixed genus. As shown in |H1 EJ, one can in fact use the orthogonal 
polynomials to provide closed expressions for F g (t) in the one-cut case. We will now 
explain how to do this in some detail. 

The object we want to compute is 

J- — F — Fq = log Z — log Zq. (2.105) 
If we write the usual series T = ^2 g>0 Fg9t 9 ~ 2 -, we have 

& = ^dog Z- log Z G) ^ log | ^ £(1 - 1) log 'jUjl. (2 , 06 ) 

The planar contribution to the free energy JF (t) is obtained from (|2.106|) by taking 
N —>■ oo. In this limit, the variable 

^ N 

becomes a continuous variable, < £ < 1, in such a way that 

N 



1 r 1 

_£/(*/JV)- / #/(0 
k=l ^° 



as A" goes to infinity. Let us assume as well that rk(N) has the following asymptotic 
expansion as A" — > oo: 

oo 

r k (N) = J2N- 2s R2s(0- (2-107) 

We then find 

F Q {t) = -h 2 log t + t 2 £ de(l - log ^Q. (2.108) 
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This provides a closed expression for the planar free energy in terms of the large N 
limit of the recursion coefficients rj,. 

It is interesting to see how to recover the density of states p(A) in the saddle-point 
approximation from orthogonal polynomials. Let us first try to evaluate (j2.1U4|) in 
the planar approximation, following jH]. A simple argument based on the recursion 
relations indicates that, at large N, 

Using now the integral representation 

l\ f 1 dy (2yY 



we find 



i(TVA/')= f M f i JL^ =(2 yRl'\ 

N Jo J-l TT J\ - y 1 



(0 - mr 

T 

where we have denoted by s(£) the limit as N — > oo of the recursion coefficients s&(iV) 
which appear in f)2.92j) . Since the above average can be also computed by (|2.52jl . by 
comparing we find 

p(A) = f\ f ^^=L=5{\-(2yRl /2 (0-s(0)), 
Jo J -i 71 V 1 ~~ V 

or, more explicitly, 

f*<mm-(\ + «m (2 . 110) 

Here, 9 is the step function. It also follows from this equation that p(A) is supported 
on the interval [b(t),a(t)], where 

b(t) = -2y/Ro(l) - s(l), a(t) = 2y/Ro(l) - s(l). (2.111) 

Example. In the Gaussian matrix model i?o(0 — t£, and s(£) = 0. We then find 
that the density of eigenvalues is supported in the interval [—2y/i, 2y/t] and it is given 
by 

^ ' Wo - A 2 2nt 

which reproduces of course Wigner's semicircle law. 
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As shown in jSJ , orthogonal polynomials can be used as well to obtain the higher 
genus free energies T g . The key ingredient to do that is simply the Euler-MacLaurin 
formula, which reads 

k=i J0 P =i v y >' 

(2.112) 

and should be regarded as an asymptotic expansion for N large which gives a way to 
compute systematically 1/N corrections. We can then use it to calculate (|2.1(J6|) at all 
orders in 1/N, where 

/(fc/iV ) = ( 1 _y iog ^W, (2 . 113 ) 

and we use the fact that has an expansion of the form ([2. 107)1 . In this way, we find 
for example that 



R (0 12 d£ 



e Jo 



and so on. We will use this formulation in section 4 to compute J-" g (t) in the matrix 
model that describes Chern-Simons theory on S 3 . 

It is clear from the above analysis that matrix models can be solved with the method 
of orthogonal polynomials, in the sense that once we know the precise form of the 
coefficients in the recursion relation we can compute all quantities in an 1/N expansion. 
Since the recursion relation is only known exactly in a few cases, we need methods 
to determine its coefficients for general potentials W(M). In the case of polynomial 
potentials, of the form 

W(M) = J^TrM p , 
p>o P 

there are well-known techniques to obtain explicit results |l] , see |21[ |2*2~] for reviews. 
We start by rewriting the recursion relation ()2.92|) as 

n+l 

APn(A) = 2j B nm p m , 
m=0 

where B is a matrix. The identities 

r n J dXe~^ ww W'(X)p n (X)p n ^(X) = nh n9s , 

J dX^-( Pn e-^ WW Pn ) = (2.114) 
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lead to the matrix equations 



is- 



{W\B)) nn _ x = ng s 
(W'(B)) nn = 0. (2.115) 

These equations are enough to determine the recursion coefficients. Consider for ex- 
ample a quartic potential 

W(A) = ^A 2 + ^A 4 . 
2 4 

Since this potential is even, it is easy to see that the first equation in (|2.115j) is auto- 
matically satisfied, while the second equation leads to 

r n {92 + 9i{r n + r n -i + r n+1 ) } = ng s 

which at large N reads 

Ro(92 + 3g 4 R ) = £t. 
In general, for an even potential of the form 

5Vf2 



one finds 



W (\) = J2 ^E±L\ 2p+2 (2.116) 
P >o ^ + 2 

^ = T,92 P J 2P+1 )r P o + \0, (2.H7) 

P >o \ P / 

which determines Rq as a function of ^. The above equation is sometimes called - 
especially in the context of double-scaled matrix models- the string equation, and by 
setting £ = 1 we find an explicit equation for the endpoints of the cut in the density of 
eigenvalues as a function of the coupling constants and t. 

Exercise. Verify, using saddle-point techniques, that the string equation correctly 
determines the endpoints of the cut. Compute i?o(£) for the quartic and the cubic 
matrix model, and use it to obtain J-o(t) (for the quartic potential, the solution is 
worked out in detail in 0). 



3 Type B topological strings and matrix models 
3.1 The topological B model 

The topological B model was introduced in ES| and can be constructed by twisting 
the Af = 2 superconformal sigma model in two dimensions. There are in fact two 
different twists, called the A and the B twist in (HUES!, and in these lectures we will 
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focus on the second one. A detailed review of topological sigma models and topological 
strings can be found in jSHl- 

The topological B model is a theory of maps from a Riemann surface T, g to a Calabi- 
Yau manifold X of complex dimension d. The Calabi-Yau condition arises in order 
to cancel an anomaly that appears after twisting (see for example Chapter 3 of [S2] 
for a detailed analysis of this issue). Indices for the real tangent bundle of X will 
be denoted by i — 1, ■ ■ ■ , 2d, while holomorphic and antiholomorphic indices will be 
denoted respectively by I , I — 1, ■ ■ ■ , d. The holomorphic tangent bundle will be simply 
denoted by TX, while the antiholomorphic tangent bundle will be denoted by TX. One 
of the most important properties of Calabi-Yau manifolds (which can actually be taken 
as their defining feature) is that they have a holomorphic, nonvanishing section Q of the 
canonical bundle Kx = Q 3,0 (X). Since the section is nowhere vanishing, the canonical 
line bundle is trivial and ci(Kx) = 0. We will always consider examples with complex 
dimension d = 3. 

The field content of the topological B model is the following. First, since it is a 
nonlinear sigma model, we have a map x : £ 9 — > X, which is a scalar, commuting field. 
Besides the field x, we have two sets of Grassmann fields t] 1 \9 I G x*(TX), which are 
scalars on E 9 , and a Grassmannian one- form on E g , p J a , with values in x*(TX). We also 
have commuting auxiliary fields F 1 , F 1 (we will follow here the off-shell formulation of 
jini]50j). The action for the theory is: 

C = t! d 2 z\G a (d z x I d- z x 7 + d- z x I d z x 7 ) - pKGjjD.r] 7 + D- Z 9 T ) 

-pi {G n D z rf - D z 9j) - R I jTK i 1 T p J zP ^e I - GjjF'F 7 ] , (3.1) 

In this action, we have picked local coordinates z, z on S 3 , and d 2 z is the measure 
—idzAdz. t is a parameter that plays the role of 1/h, the field 0j is given by 6>/ = G I -j9 J , 
and the covariant derivative D a acts on sections ip l of the tangent bundle as 

Datf = d a f + d a x>ry k . (3.2) 

The theory also has a BRST, or topological, charge Q which acts on the fields according 
to 

= T] 1 , 
1 G IJ FJ i 

-- D zP J -_D.pl + R I jTK r l L p J z pf, 

-. -T I —n J F K 
jk I ' 

The action of Q explicitly depends on the splitting between holomorphic and anti- 
holomorphic coordinates on X, in other words, it depends explicitly on the choice of 



[Q,x T ] 


= o, 


[Q,x ! ] 


{Q,v r } 


= o, 


{Q,0i} 


{Q,pI} 


= d z x J , 




{Q,pI} 


= dz-x 1 , 
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complex structure on X. It is easy to show that Q 2 = 0, and that the action of the 
model is Q-exact: 

C = {Q,V} (3.3) 
where V (sometimes called the gauge fermion) is given by 



V = t [ d 2 z [Gjj(pid s x J + pld z x J ) - F'Oj] . 



(3.4) 



Finally, we also have a U(l) ghost number symmetry, in which x, rj, 9 and p have 
ghost numbers 0, 1, 1, and —1, respectively. The Grassmannian charge Q then has 
ghost number 1. Notice that, if we interpret rj 1 as a basis for antiholomorphic dif- 
ferential forms on X, the action of Q on x 1 , x 1 may be interpreted as the Dolbeault 
antiholomorphic differential d. 

It follows from ()3.3j) that the energy-momentum tensor of this theory is given by 

T af3 = {Q,b af3 }, (3.5) 

where b a p = 5V/5g al3 and has ghost number —1. The fact that the energy- momentum 
tensor is Q-exact means that the theory is topological, in the sense that the partition 
function does not depend on the background two-dimensional metric. This is easily 
proved: the partition function is given by 

£>0e~ £ , (3.6) 

where denotes the set of fields of the theory, and we compute it in the background of 
a two-dimensional metric g a p on the Riemann surface. Since T a @ = 5C/5g a ^, we find 
that 

<)Z <{Q,M>, (3.7) 



5g a P 

where the bracket denotes an unnormalized vacuum expectation value. Since Q is a 
symmetry of the theory, the above vacuum expectation value vanishes, and we find 
that Z is metric- independent, at least formally. 

The Q-exactness of the action itself also has an important consequence: the same 
argument that we used above implies that the partition function of the theory is inde- 
pendent of t. Now, since t plays the role of 1/h, the limit of t large corresponds to the 
semiclassical approximation. Since the theory does not depend on t, the semiclassical 
approximation is exact. The classical configurations for the above action are constant 
maps x : S 9 — > X. Therefore, it follows that path integrals of the above theory reduce 
to integrals over X [73*] . 
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What are the operators to consider in this theory? Since the most interesting aspect 
of this model is the independence w.r.t. to the two-dimensional metric, we want to 
look for operators whose correlation functions satisfy this condition. It is easy to see 
that the operators in the cohomology of Q do the job: topological invariance requires 
them to be Q-closed, and on the other hand they cannot be Q-exact, since otherwise 
their correlation functions would vanish. One can also check that the Q-cohomology 
is given by operators of the form 

O^^-V 1 '--^---^, (3-8) 

where 

= d'-^dx 11 A • ■ ■ A dx 1 * — -r A • • ■ A — -r (3.9) 

is an element of H^(X, A q TX) . Therefore, the Q-cohomology is in one-to-one corre- 
spondence with the twisted Dolbeault cohomology of the target manifold X. We can 
then consider correlation functions of the form 

(U°^- ( 3 - 10 ) 

a 

This correlation function vanishes unless the following selection rule is satisfied 

1>« = X> = <!-</), (3.11) 

a a 

where g is the genus of the Riemann surface. This selection rule comes from a 
U{1)l x U{1)r anomalous global current. Due to the arguments presented above, 
this correlation function can be computed in the semiclassical limit, where the path 
integral reduces to an integration over the target X. The product of operators in 
corresponds to a form in H^(X, A d TX) . To integrate such a form over X we crucially 
need the Calabi-Yau condition. This arises as follows. In a Calabi-Yau manifold we 
have an invertible map 

Q°' p (A q TX) — > Q d - g ' p (X) 

where the (d, 0)-form Q is used to contract the indices. Since Q is holomorphic, this 
descends to the <9-cohomology. It then follows that an element in H^(X, A d TX) maps 
to an element in H^(X). After further multiplication by Q, one can then integrate 
a (d, (i)-form over X. This is the prescription to compute correlation functions like 
()3.10|) . A simple and important example of this procedure is the case of a Calabi- 
Yau threefold, d = 3, and operators associated to forms in H±(X,TX), or by using 
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([3. 12)) . to forms in (X). These operators are important since they correspond to 
infinitesimal deformations of the complex structure of X. The selection rule (|3.11|) says 
that we have to integrate three of these operators, and the correlation function reads 
in this case 

(o^o^ 3 ) = I (&)£(&) J (0 3 )| n Ilhh dz^dz l2 dz^ a a (3.13) 

It turns out that the full information of the correlators ()3.13|1 at genus zero can 
be encoded in a single function called the prepotential. We will quickly review here 
some of the basic results of special geometry and the theory of the prepotential for 
the topological B model, and we refer the reader to [T71 EH] for more details. The 
correlation functions in the B model, like for example ()3.13)1 . depend on a choice of 
complex structure, as we have already emphasized. The different complex structures 
form a moduli space Ai of dimension h 2 ' 1 . A convenient parametrization of M. is 
the following. Choose first a symplectic basis for H 3 (X), denoted by (A a ,B a ), with 
a = 0, 1, ■ • • , h 2 ' 1 , and such that A a n B b = 5 b a . We then define the periods of the 
Calabi-Yau manifold as 

z a = [ Q, T a = [ Q, a = (),••• ,/i 2 ' 1 . (3.14) 

J A a J B a 

Of course, the symplectic group Sp(2h 2 ' 1 + 2,R) acts on the vector (^ a ,jF a ). A basic 
result of the theory of deformation of complex structures says that the z a are (locally) 
complex projective coordinates for Ai. Inhomogeneous coordinates can be introduced 
in a local patch where one of the projective coordinates, say zq, is different from zero, 
and taking 

t a = ^, a = l,.-. ,h 2 >\ (3.15) 
Zq 

The coordinates z a are called special projective coordinates, and since they parametrize 
M. we deduce that the other set of periods must depend on them, i.e. T a = T a {z). 
Using the periods ()3.14|) we can define the prepotential T{z) by the equation 



The prepotential satisfies 



T= X -Y^z a T a . (3.16) 

a=0 

*•(,) = g (3.17) 

and turns out to be a homogeneous function of degree two in the z a . Therefore, one 
can rescale it in order to obtain a function of the inhomogeneous coordinates t a : 

F (t a ) = \j r (Za). (3-18) 
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The fact that the coordinates z a are projective is related to the freedom in normaliz- 
ing the three-form Q. In order to obtain expressions in terms of the inhomogeneous 
coordinates t a , we simply have to rescale fl — > — fl, and the periods (z a , become 

(M.,^-E*.f,f). (3.19) 

a=l 

One of the key results in special geometry is that the correlation functions ()3.13j) can 
be computed in terms of the prepotential Fo(t). Given a deformation of the complex 
structure parametrized by t a , the corresponding tangent vector d/dt a is associated to a 
differential form of type (2, 1). This form leads to an operator O a , and the three-point 
functions involving these operators turn out to be given by 

(O a O b O c ) = - \ ° . (3.20) 
x ' dt a dt b dt c y ' 

The prepotential F (t) encodes the relevant information about the B model on the 
sphere, and it has an important physical meaning, since it gives the four-dimensional 
supergravity prepotential of type IIB string theory compactified on X (and determines 
the leading part of the vector multiplet effective action). 

In order to obtain interesting quantities at higher genus one has to couple the topo- 
logical B model to two-dimensional gravity, using the fact that the structure of the 
twisted theory is very close to that of the bosonic string [301 EH EZj • In the bosonic 
string, there is a nilpotent BRST operator, Qbrst, and the energy-momentum tensor 
turns out to be a QBRST-commutator: T(z) = {Qbrst, K z )}- m addition, there is a 
ghost number with anomaly 3x(S 9 ) = 6 — 6g, in such a way that Qbrst and b(z) have 
ghost number 1 and —1, respectively. This is precisely the same structure that we 
found in ()3.5j) . and the composite field b a/ 3 plays the role of an antighost. Therefore, 
one can just follow the prescription of coupling to gravity for the bosonic string and 
define a genus g > 1 free energy as follows: 

F v= /_ (II ( & '^))' ( 3 - 21 ) 

JM 3 k=l 

where 

(6,^)=/ dMM^/ + M74)/), (3-22) 
Jt, 9 

and fik are the usual Beltrami differentials. The vacuum expectation value in (|3.21j) 
refers to the path integral over the fields of the topological B model, and gives a 
differential form on the moduli space of Riemann surfaces of genus g, M g , which is 
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then integrated over. The free energies F g of the B model coupled to gravity for g > 1 
are also related to variation of complex structures. A target space description of this 
theory, called Kodaira-Spencer theory of gravity, was found in 7], and can be used to 
determine recursively the F g in terms of special geometry data. 

3.2 The open type B model and its string field theory descrip- 
tion 

The topological B model can be formulated as well for open strings, i.e., when the 
worldsheet is an open Riemann surface with boundaries jZH EH]- In order to 
construct the open string version we need boundary conditions (b.c.) for the fields. It 
turns out that the appropriate b.c. for the B model are Dirichlet along holomorphic 
cycles of X, S, and Neumann in the remaining directions. Moreover, one can add Chan- 
Paton factors to the model, and this is implemented by considering a U(N) holomorphic 
bundle over the holomorphic cycle S. The resulting theory can then be interpreted as 
a topological B model in the presence of N topological D-branes wrapping S. Since we 
will be interested in finding a spacetime description of the open topological B model, 
we can consider the case in which the branes fill spacetime (the original case considered 
in (71]) and deduce the spacetime action for lower dimensional branes by dimensional 
reduction. In the spacetime filling case, when S = X, the boundary conditions for the 
fields are 9 = along dUg^ and that the pullback to <9£ 9) h of *p vanishes (where * is 
the Hodge operator). 

The open topological B model can also be coupled to gravity following the same 
procedure that is used in the closed case, and one obtains in this way the open type B 
topological string propagating along the Calabi-Yau manifold X. We are now interested 
in providing a description of this model when the N branes are spacetime filling. As 
shown by Witten in |73| . the most efficient way to do that is to use the cubic string 
field theory introduced in [69J. 

In bosonic open string field theory we consider the worldsheet of the string to be an 
infinite strip parameterized by a spatial coordinate < a < tt and a time coordinate 
— oo < t < oo, and we pick the flat metric ds 2 = da 2 + dr 2 . We then consider 
maps x : I —> X, with / = [0, tt] and X the target of the string. The string field 
is a functional of open string configurations \l/[a;(cr)], of ghost number one (the string 
functional depends as well on the ghost fields, but we do not indicate this dependence 
explicitly). In Witten defines two operations on the space of string functionals. 
The first one is the integration, which is defined formally by folding the string around 
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its midpoint and gluing the two halves: 




(3.23) 



The integration has ghost number —3, which is the ghost number of the vacuum. This 
corresponds to the usual fact that in open string theory on the disk one has to soak 
up three zero modes. One also defines an associative, noncommutative star product * 
of string functionals through the following equation: 




where xn+i = x\. The star product simply glues the strings together by folding them 
around their midpoints, and gluing the first half of one with the second half of the 
following (see for example the review jHHj for more details), and it doesn't change the 
ghost number. In terms of these geometric operations, the string field action is given 
by 

S = j J Qtf * Qbrst^ + itf * * J (3.25) 

where g s is the string coupling constant. Notice that the integrand has ghost number 
3, while the integration has ghost number —3, so that the action (J3.25|) has zero ghost 
number. If we add Chan-Paton factors, the string field is promoted to a U(N) matrix 
of string fields, and the integration in (J3.25|) includes a trace Tr. The action (J3.25j) 
has all the information about the spacetime dynamics of open bosonic strings, with or 
without D-branes. In particular, one can derive the Born-Infeld action describing the 
dynamics of D-branes from the above action [S3] . 

We will not need all the technology of string field theory in order to understand 
open topological strings. The only piece of relevant information is the following: the 
string functional is a function of the zero mode of the string (which corresponds to 
the position of the string midpoint), and of the higher oscillators. If we decouple 
all the oscillators, the string functional becomes an ordinary function of spacetime, 
the * product becomes the usual product of functions, and the integral is the usual 
integration of functions. The decoupling of the oscillators is in fact the point-like limit 
of string theory. As we will see, this is the relevant limit for topological open type B 
strings on X. 

We can now exploit again the analogy between open topological strings and the open 
bosonic string that we used to define the coupling of the topological B model to gravity 
(i.e., that both have a nilpotent BRST operator and an energy-momentum tensor that 
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is QBRST-exact). Since both theories have a similar structure, the spacetime dynamics 
of open topological type B strings is governed as well by (|3.25|) . where Qbrst is given 
in this case by the topological charge defined in ()3.3|) . and where the star product and 
the integration operation are as in the bosonic string. The construction of the cubic 
string field theory also requires the existence of a ghost number symmetry, which is 
also present in the topological sigma model in the form of a U(1)r symmetry, as we 
discussed in 3.1. It is convenient to consider the U(1)r charge of the superconformal 
algebra in the Ramond sector, which is shifted by —d/2 with respect to the assignment 
presented in 3.1 (here, d is the dimension of the target). When d = 3 this corresponds 
to the normalization used in 69], in which the ghost vacuum of the be system is assigned 
the ghost number —1/2. 

In order to provide the string field theory description of open topological type B 
strings on X, we have to determine the precise content of the string field, the * algebra 
and the integration of string functionals for this particular model. As in the conven- 
tional string field theory of the bosonic string, we have to consider the Hamiltonian 
description of topological open strings. We then take S to be an infinite strip and con- 
sider maps x : I — > X, with I = [0, tt]. The Hilbert space is made up out of functionals 
ty[x(a), ■ ■ ■}, where x is a map from the interval as we have just described, and the • • • 
refer to the Grassmann fields (which play here the role of ghost fields). Notice that, 
since p T ZI are canonically conjugate to r], 9, we can choose our functional to depend 
only on rj, 9. It is easy to see that the Hamiltonian has the form 



We then see that string functionals with dx 1 /da ^ cannot contribute: as we saw in 
the previous subsection, the physics is t-independent, therefore we can take t — > oo. 
In this limit the functional gets infinitely massive and decouples from the spectrum, 
unless dx 1 /da = 0. Therefore, the map x : I —>■ X has to be constant and in particular 
it must be a point in X. A similar analysis holds for the Grassmann fields as well. 
Since 9 = at the boundary, it follows that string functionals are functions of the 
commuting zero modes x % and rj 1 , and can be written as 



These functionals can be interpreted as a sum of (0,p)-forms on X. If we have N D- 
branes wrapping X, these forms will be valued in End(.E') (where E is a holomorphic 
U (N) bundle). The Q symmetry acts as on these functionals as the Dolbeault operator 
d with values in End(E). Notice that a differential form of degree p will have ghost 
number p. 




(3.26) 




(3.27) 



p>i 
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We are now ready to write the string field action for topological open type B strings 
on X with N spacetime filling branes. We have seen that the relevant string functionals 
are of the form ([3.27)1 . Since in string field theory the string field has ghost number 
one, we must have 

= rfA T (x), (3.28) 

where Aj(x) is a (0, l)-form taking values in the endomorphisms of some holomorphic 
vector bundle E. In other words, the string field is just the (0, 1) piece of a gauge 
connection on E. Since the string field only depends on commuting and anticommuting 
zero modes, the star product becomes the wedge products of forms in fi(°' p )(End(.E)), 
and the integration of string functionals becomes ordinary integration of forms on X 
wedged by Q. We then have the following dictionary: 

* -> A, Qbrst -> d 

(3.29) 

The string field action ([3.25)1 is then given by 

S = — f A Tr ( A A dA + -A A A A A 
29s Jx V 3 

This is the so-called holomorphic Chern-Simons action. It is a rather peculiar quan- 
tum field theory in six dimensions, but as we will see, when we consider D-branes of 
lower dimension, we will be able to obtain from ()3.30|) more conventional theories by 
dimensional reduction. 

3.3 Topological strings and matrix models 

We have seen that the spacetime description of the open B model with spacetime 
filling branes reduces to a six-dimensional theory ([3.30)1 . We will see now that, in some 
circumstances, this theory simplifies drastically and reduces to a matrix model. 

In order to simplify the spacetime description one should study simple Calabi-Yau 
manifolds. The simplest example of a local Calabi-Yau threefold is a Riemann surface 
together with an appropriate bundle over it. The motivation for considering this kind 
of models is the following. Consider a Riemann surface S g holomorphically embedded 
inside a Calabi-Yau threefold X, and let us consider the holomorphic tangent bundle 
of X restricted to S s . We then have 

TX\z g = TS 9 © N^ g (3.31) 



(3.30) 
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where N% is a holomorphic rank two complex vector bundle over E fl , called the normal 
bundle of S g , and the CY condition Ci(X) = gives 

c x (Nv g ) = 2g-2. (3.32) 

The Calabi-Yau X "near S 5 " looks precisely like the total space of the bundle 

N -> S g (3.33) 

where A/" is regarded here as a bundle over S 5 satisfying ()3.32|) . The space ()3.33|) is an 
example of a local Calabi-Yau threefold, and it is noncompact. 

When g = and S g = P 1 it is possible to be more precise about the bundle N. A 
theorem due to Grothendieck says that any holomorphic bundle over P 1 splits into a 
direct sum of line bundles (for a proof, see for example j^Sl, pp. 516-7). Line bundles 
over P 1 are all of the form 0(n), where n e Z. The bundle 0(n) can be easily 
described in terms of two charts on P 1 : the north pole chart, with coordinates z, $ for 
the base and the fiber, respectively, and the south pole chart, with coordinates z', 
The change of coordinates is given by 

z ' = \j z , $' = Z - n <$>. (3.34) 

We also have that c\(0(n)) = n. We then find that local Calabi-Yau manifolds that 
are made out of a two-sphere together with a bundle over it are all of the form 

0(-a)@0(a-2) P 1 , (3.35) 

since the degree of the bundles have to sum up to —2 due to ()3.32|) . 

Let us now consider the string field theory of type B open topological strings on 
the Calabi-Yau manifold ()3.35|) . We will consider a situation where we have Dirichlet 
boundary conditions associated to P 1 , in other words, there are A" topological D-branes 
wrapping P 1 . Since the normal directions to the D-brane worldvolume are noncompact, 
the spacetime description can be obtained by considering the dimensional reduction 
of the original string field theory action ()3.3()j) . As usual in D-brane physics, the 
gauge potential A splits into a gauge potential on the worldvolume of the brane and 
Higgs fields describing the motion along the noncompact, transverse directions. In a 
nontrivial geometric situation like the one here, the Higgs fields are sections of the 
normal bundle. We then get three different fields: 

A $o, $i, (3.36) 

where A is a U(N) (0, 1) gauge potential on P 1 , $ is a section of 0(— a), and $i is a 
section of Oia — 2). Both fields, $o and $i, take values in the adjoint representation 



38 



of U(N). It is easy to see that the action ()3.30|) becomes 

S= - I Tz(® D A ®x), (3.37) 

where Da = d+ [A, ■] is the antiholomorphic covariant derivate. Notice that this theory 
is essentially a gauged system, since $ , 3>i are quasiprimary conformal fields of 
dimensions a/2, 1 — a/2, respectively. 

We will now consider a more complicated geometry. We start with the Calabi-Yau 
manifold ()3.35|) with a = 0, i.e. 

£>(0)©£>(-2) ^P 1 . (3.38) 

In this case, $o is a scalar field on P , while $i is a (1, 0) form (since K F i = 0(—2)). 
If we cover P 1 with two patches with local coordinates z,z' related by z' = 1/z, the 
fields in the two different patches, $0) ^i? an d $' , will be related by 

$' = $ , = ^ 2 $i. (3.39) 

We can regard this geometry as a family of P x s located at $' x = (the zero section of 
the nontrivial line bundle 0{— 2)) parametrized by $o — = x £ C. The idea is to 
obtain a geometry where we get n isolated P x s at fixed positions of x. To do that, we 
introduce an arbitrary polynomial of degree n + 1 on $ > ^(^o)> an d we modify the 
gluing rules above as follows [To] : 

z > = 1/ Zj $' Q = $ , $; = z 2 ^ + W'(<$> )z. (3.40) 

Before, the P 1 was in a family parameterized by $o £ C. Now, we see that there are 
n isolated P x s located at fixed positions of $ given by W'(Q>o) = 0, since this is the 
only way to have $i = = 0. 

The geometry obtained by imposing the gluing rules (|3.4U|) can be interpreted in yet 
another way. Call $o = % and define the coordinates 

M = 2$;, u = 2$i, y = i{2z'& 1 -W'{x)). (3.41) 

The last equation in (j3.40j) can now be written as 

uv + y 2 + W\x) 2 = 0. (3.42) 

This is a singular geometry, since there are singularities along the line u = v = y = 
for every such that W'(x*) = 0. For example, if W'(x) = x, (13. 42 J) becomes, after 
writing u,v u — iv,u + iv 

v 2 + v 2 + x 2 + y 2 = 0. (3.43) 
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This Calabi-Yau manifold is called the conifold, and it is singular at the origin. For ar- 
bitrary polynomials W(x), the equation ([3.42)1 describes more general, singular Calabi- 
Yau manifolds. Notice that locally, around the singular points u — v — y — 0, x = x*, 
the geometry described by (j3.42j) looks like a conifold (whenever W"(x*) = 0). The 
manifold described by ()3.40j) is obtained after blowing up the singularities in ()3.42|) . 
i.e. we modify the geometry by "inflating" a two-sphere IP 1 at each singularity. This 
process is called resolution of singularities in algebraic geometry, and for this reason 
we will call the manifold specified by (j3.40j) the resolved manifold X rcs . 

We can now consider the dynamics of open type B topological strings on X res . We 
will consider a situation in which we have in total N D-branes in such a way that Ni 
D-branes are wrapped around the i-th IP 1 , with i — 1, • • ■ , n. As before, we have three 
fields in the adjoint representation of U(N), $0; an d the gauge connection A. The 
action describing the dynamics of the D-branes turns out to be given by 

S = — [ Tr ( $iD A $o + wW($ ) J (3.44) 
9s Jpi V / 

where a; is a Kahler form on P 1 with unit volume. This action was derived in Ol^Zj. A 
quick way to see that the modification of the gluing rules due to adding the polynomial 
W'($o) leads to the extra term in (|3.44|) is to use standard techniques in CFT |27| . 
The fields <3?o ; are canonically conjugate and on the conformal plane they satisfy the 
OPE 

$ (z)$i(«j) — . (3.45) 

z — w 

Let us now regard the geometry described in ()3.40|) as two disks (or conformal planes) 
glued through a cylinder. Since we are in the cylinder, we can absorb the factors of z 
in the last equation of ()3.40|) . The operator that implements the transformation of $ 
is 

U = exp — / Tr W(Q (z)) dz, (3.46) 
9s J 

since from 1)3.45)1 it is easy to obtain 

$; = U^U- 1 . (3.47) 

We can also write 

f/ = exp— / TW($ (z)V (3.48) 
9s Jv 1 

where uj is localized to a band around the equator of P 1 (as we will see immediately, 
the details of uj are unimportant, as long as it integrates to 1 on the two-sphere). 
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One easy check of the above action is that the equations of motion lead to the 
geometric picture of D-branes wrapping n holomorphic P x s in the geometry. The 
gauge connection is just a Lagrange multiplier enforcing the condition 

[$o,$i] = 0, (3.49) 

therefore we can diagonalize $o and $1 simultaneously. The equation of motion for $o 
is simply 

M Q = 0, (3.50) 

and since we are on P , we have that <3>o is a constant, diagonal matrix. Finally, the 
equation of motion for $! is 

M x = W'(<$> )u, (3.51) 

and for nonsingular $i configurations both sides of the equation must vanish simulta- 
neously, as we can see by integrating both sides of the equation over P 1 . Therefore, 
$i = and the constant eigenvalues of $o satisfy 

W'(%) = (3.52) 

i.e. they must be located at the critical points of W(x). In general, we will have 
Ni eigenvalues of $o & t the i-th critical point, % = 1, • • • ,n, and this is precisely the 
D-brane configuration we are considering. 

What happens in the quantum theory? In order to analyze it, we will use the 
approach developed in ^H] for the analysis of two-dimensional gauge theories 2 . First 
of all, we choose the maximally Abelian gauge for <3> , i.e. we write 

$ = ^ + $*, (3.53) 

where $q is the projection on the Cartan subalgebra t, and $g is the projection on the 
complementary part k. The maximally Abelian gauge is defined by the condition 

<5>o=0 (3.54) 

which means that the nondiagonal entries of $o are gauge-fixed to be zero. This is in 
fact the same gauge that we used before to write the matrix model in the eigenvalue 
basis. After fixing the gauge the usual Faddeev-Popov techniques lead to a ghost 
functional determinant given by 

^Det k (ad($*)) n0( pi ) (3.55) 

2 I'm grateful to George Thompson for very useful remarks on this derivation. 
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where the subscript k means that the operator $q acts on the space k, and the normal- 
ization factor 1/N\ is the inverse of the order of the residual symmetry group, namely 
the Weyl group which permutes the N entries of $q. The integrand of (|3.44j) reads, 
after gauge fixing, 

Tr (^d^l + W^l) J +2^A a ^ a a{$l), (3.56) 

\ 'a 

where a are roots, E a is a basis of k, and we have expanded $1 = S Q ^\E a as well 
as A k . We can now integrate out the A a to obtain 

, 77-rr, TT ( 3 - 57 ) 

Det k (ad($ t )) nl ,o (pl) 11 V 1} V } 

Here we have used the functional generalization of the standard formula 5 (ax) = 
\a\~ 1 8(x). We can now trivially integrate over The inverse determinant in ()3.57j) 
combines with ()3.55|) to produce 

Det k (ad($ t )) g o (P i ) 
Det k (ad($J)) H i,o(pi) 

where (as usual) nonzero modes cancel (since they are paired by d) and one ends with 
the determinants evaluated at the cohomologies. Similarly, integrating out in (|3.56|) 
leads to <9$q = 0, therefore $q must be constant. The quotient of determinants is easy 
to evaluate in this case, and one finds 

h°(P 1 )-h 1 -°(p 1 ) 



i<j 



(3.59) 



where Aj are the constant eigenvalues of $q. Since ^(P 1 ) = 1, /i 1,0 (P 1 ) = 0, we just 
get the square of the Vandermonde determinant and the partition function reads: 

i r N 

' J i=l i<j 

In principle, as explained in ^U], one has to include a sum over nontrivial topological 
sectors of the Abelian gauge field A* in order to implement the gauge fixing ()3.54j) 
correctly. Fortunately, in this case the gauge- fixed action does not depend on A 1 , and 
the inclusion of topological sectors is irrelevant. The expression (|3.6U|) is (up to a factor 
(271*)^) the gauge-fixed version of the matrix model 



jv^e-fs™^ (3.61) 



vo\(U(N)) 

We have then derived a surprising result due to Dijkgraaf and Vafa ^T7\: the string 
field theory action for open topological B strings on the Calabi-Yau manifold described 
by (|3.4U|) is a matrix model with potential W(&). 
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3.4 Open string amplitudes and multicut solutions 



The total free energy F(Ni,g s ) of topological B strings on the Calabi-Yau (|3.40|) in 
the background of N = ^iVj branes wrapped around n P^s is of the form (|1.3j) . 
and as we have just seen it is given by the free energy of the matrix model ()3.61|) . 
In particular, the coefficients Fg^,- ,h„ can be computed perturbatively in the matrix 
model. We have to be careful however to specify the classical vacua around which we are 
doing perturbation theory. Remember from the analysis of the matrix model that the 
classical solution which describes the brane configuration is characterized by having iV, 
eigenvalues of the matrix located at the i-th critical point of the potential W(x). In the 
saddle-point approximation, this means that we have to consider a multicut solution, 
with eigenvalues "condensed" around all the extrema of the potential. Therefore, in 
contrast to the multicut solution discussed in 2.2, we have that (1) all critical points 
of W(x) have to be considered, and not only the minima, and (2) the number of 
eigenvalues in each cut is not determined dynamically as in ([2.83)1 . but it is rather 
fixed to be Ni in the i-th cut. In other words, the integral of the density of eigenvalues 
p(X) along each cut equals a fixed filling fraction z/j = Ni/N: 



where N = Yli=i -^i * s the total number of eigenvalues. Let us introduce the partial 't 
Hooft couplings 



where Ai is the closed cycle of the hyperelliptic curve ()2.7()|) which surrounds the cut Cj. 
Assuming for simplicity that all the £j are different from zero, and taking into account 
that J2i U = we see that (|3.fi4j) gives n — 1 independent conditions, where n is the 
number of critical points of W(x). These conditions, together with (|2.81jh determine 
the positions of the endpoints Xi as functions of the £j and the coupling constants in 
W(x). It is clear that the solution obtained in this way is not an equilibrium solution 
of the matrix model, since cuts can be centered around local maxima and different 
cuts will have different values of the effective potential. This is not surprising, since 
we are not considering the matrix model as a quantum mechanical system per se, but 
as an effective description of the original brane system. The different choices of filling 
fractions correspond to different choices of classical vacua for the brane system. 




(3.62) 



U = g s Ni = tvi. 



(3.63) 



Taking into account ()2.58|) and ()2.68|) . we can write ()3.62|) as 




(3.64) 
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A subtle issue concerning the above matrix model is the following. The matrix field 
$ in (|3.61|) comes from the B model field $o> which is a holomorphic field. Therefore, 
the matrix integral ()3.60|) should be understood as a contour integral, and in order to 
define the theory a choice of contour should be made. This can be done in perturbation 
theory, by choosing for example a contour that leads to the usual results for Gaussian 
integration, and therefore at this level the matrix model is not different from the 
usual Hermitian matrix model [23 ESI- I n some cases, however, regarding (j3.61|) as a 
holomorphic matrix model can be clarifying, see [5TJ for an exhaustive discussion. 

The above description of the multicut solution refers to the saddle-point approx- 
imation. What is the meaning of the multicut solutions from the point of view of 
perturbation theory? To address this issue, let us consider for simplicity the case of 
the cubic potential: 

— W($) = — Tr$ 2 + -— Tr$ 3 . (3.65) 
9s ^9s 3 g s 

This potential has two critical points, a\ = and a 2 = — 1/(3. The most general 

multicut solution will have two cuts. There will Ni eigenvalues sitting at $ = 0, and 

N2 eigenvalues sitting at $ = — 1/(3. The partition function Z of the matrix model is: 



i=l 



where A(A) = rL<j(^« — ^j) * s the Vandermonde determinant. We can now expand 
the integrand around the vacuum with Aj = for i = 1, . . . , iVi and A, = — -5 for % = 
Ni + 1, . . . , N. Denoting the fluctuations by /ij and Uj, the Vandermonde determinant 
becomes 

a 2 (a)= n ta-/o a n (yh-»h) a n (^-^ + l) 2 - ^ 

Kii<i 2 <ATi Kji<j2<N 2 i<i<Ni " 

l<j<N 2 

We also expand the potential around this vacuum and get 

Notice that the propagator of the fluctuations around — 1/(3 has the 'wrong' sign, since 
we are expanding around a local maximum. The interaction between the two sets 
of eigenvalues, which is given by the last factor in ()3.67|) . can be exponentiated and 
included in the action. This generates an interaction term between the two eigenvalue 
bands 

W int = 2N,N 2 log (3 + 2 £ i/3 fc ^(-1)'* fiJf*. (3.69) 

k=l i,j p=0 
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By rewriting the partition function in terms of matrices instead of their eigenvalues, 
we can represent this model as an effective two-matrix model, involving an N\ x Ni 
matrix $ 1; and an N 2 x N 2 matrix $ 2 : 



1 



Yo\(U{N 1 )) x Vo\(U(N 2 )) 



J D^iD^2e' Wl ^ l) ~ W2 ^ 2) ' w ^ 1 ' q " 2 \ (3.70) 



where 



W 2 {<$> 2 ) 
W int ($ u $ 2 ) 



+Trf— $? 



-Trf^^--^,, 



2<? s 
k 



+ A$3 

3<? s 1 
3^ 



Tr Tr $ 2 ? 



fc=l p =o 

+ iV 2 W(a 2 ) + iV 1 W r (a 1 ) -2^!iV 2 ln/5 



k—p 



(3.71) 



Here, Tr$° = JVi, Tr$° = JV 2 , W(ai) = and W(a 2 ) = l/{6g s /3 2 ). Although the 
kinetic term for $ 2 has the 'wrong' sign, we can still make sense of the model in 
perturbation theory by using formal Gaussian integration, and this can in fact be 
justified in the framework of holomorphic matrix models [HI]. Therefore, the two- 
cut solution of the cubic matrix model can be formally represented in terms of an 
effective two-matrix model. It is now straightforward to compute the free energy F peTt = 
\og(Z((3)/Z((3 = 0)) in perturbation theory. It can be expanded as 



F pe * = -N 1 W(a 1 )-N 2 W(a 2 )-2N 1 N 2 \nP + J2J2^sP 2 ) 29 ~ 2+h F a ,h(N 1 ,N 2 ) (3.72) 

h=l g>0 

where F g ^ is a homogeneous polynomial in Ni and N 2 of degree h. One finds, up to 
fourth order in the coupling constant f3, the following result 



+ 9sP 



pert 
2 



-iVW(ai) - N 2 W{a 2 ) - 2iV 1 iV 2 ln/3 



-N? - 5N*N 2 + 5N X N 2 - -N$ 

3 1 12 2 3 2 



+ -(N 1 -N 2 ) 
o 



+ & 
+ ■■■ 



(In? - yiVfiV 2 + 59iVM - 9 -± Nl Ni + 8 -N}) + (IjV? - ^ Nl N 2 + IjVf 



(3.73) 



From this explicit perturbative computation one can read off the first few coefficients 
FgMM- 0^ course, this procedure can be generalized, and the n-cut solution can 
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be represented by an effective n matrix model with interactions among the different 
matrices that come from the expansion of the Vandermonde determinant. These in- 
teractions can be also incorporated in terms of ghost fields, as explained in [21]. This 
makes possible to compute corrections to the saddle-point approximation in perturba- 
tion theory. One can also use the multicut solution to the loop equations jH E7| with 
minor modifications to compute the genus one correction in closed form |1H | I26 | IT%]. 



3.5 Master field and geometric transition 

We have seen that the open topological string amplitudes on the Calabi-Yau manifold 
X res are computed by a multicut matrix model whose planar solution (or, equivalently, 
its master field configuration) is given by a hyperelliptic curve 

y 2 = W'(x) 2 -R(x). (3.74) 

Moreover, we also saw in ()3.64|) that the partial 't Hooft couplings can be understood 
as integrals around the A± cycles of this curve, with i = 1, ■ ■ • , n. Let us now compute 
the variation of the free energy F (ti) when we vary tj. The variation w.r.t. ti (keeping 
the tj, j i, fixed) can be obtained by computing the variation in the free energy as we 
move one eigenvalue from the cut Ci to infinity |27j . This variation is given by (minus) 
the integral of the force exerted on an eigenvalue, as we move it from the endpoint 
of the cut to infinity. The path from the endpoint of Ci to infinity, which does not 
intersect the other cuts Cj, will be denoted by Bi. Taking into account f)2.69|) . and the 
fact that y(p) has no discontinuities outside the cuts Cj, we find 



dFo 
dt 



y(x)dx. (3.75) 



Usually this integral is divergent, but can be easily regularized by taking Bi to run up 
to a cutoff point x = A, and subtracting the divergent pieces as the cutoff A goes to 
infinity. For example, for the Gaussian matrix model one has 



— -=/ dx Vx 2 -At = t(log t - 1) - 2t log A + -A 2 + C(l/A 2 ). (3.76) 
vt J2VI 2 



Therefore, the regularized integral gives t(log t — 1), which is indeed the right result. It 
is now clear that (|3.64j) and (|3.75|) look very much like the relations (|3.14j) that define 
the periods (therefore the prepotential) in special geometry. What is the interpretation 
of the appearance of special geometry? 

Recall that our starting point was a Calabi-Yau geometry obtained as a blowup of the 
singularity given in (j3.42j) . However, there is another way of smoothing out singularities 
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in algebraic geometry, which is by deforming them rather than by resolving then. For 
example, the conifold singularity given in (J3.43)) can be smoothed out by deforming 
the geometry to 

x 2 + y 2 + u 2 + v 2 = fi. (3.77) 

This is the so called deformed conifold. Geometrically, turning on /j corresponds to 
inflating a three-sphere in the geometry, since the real section of the conifold is indeed 
an S 3 . As /x — > 0, the three-sphere collapses to zero size, so we can interpret the 
singularity as arising from a collapsing three-cycle in the geometry. In the more general 
singularity 1)3.42)1 . the generic deformation requires turning on a generic polynomial of 
degree n — 1 R(x), and we get the Calabi-Yau manifold 

u 2 + v 2 + y 2 + W'(x) 2 = R(x). (3.78) 

We will call this geometry the deformed manifold X^ef- The deformation by R(x) 
introduces in fact n three-spheres in the geometry, one for each singularity (recall that 
each of the singular points in (J3.42)) is locally like the conifold). The noncompact 
Calabi-Yau manifold (J3.78)) has a holomorphic three-form: 

n = (3.79) 

2n v 

The three-spheres created by the deformation can be regarded as two-spheres fibered 
over an interval in the complex x-plane. To see this, let us consider for simplicity the 
case of the deformed conifold ([3.77)1 . with /i real. This geometry contains a three- 
sphere which is given by the restriction of ()3.77j) to real values of the variables. If 
we now consider a fixed, real value of x in the interval — JJi < x < y/Ji, we get of 
course a two-sphere of radius \J \i — x 2 . The sphere collapses at the endpoints of the 
interval, x = ±-y//I, and the total geometry of the two-sphere together with the interval 
[— y/JI, JJj\ is a three-sphere. In the more general case, the curve W'(x) 2 — R(x) has n 
cuts with endpoints £21,2^-1, i = 1, • • • , n , and the n three-spheres are S 2 fibrations 
over these cuts. 

Let us now consider closed type B topological strings propagating on X^ e {. As we 
saw in 3.1, the genus zero theory is determined by the periods of the three- form Q given 
in (|3.79j) . We then choose a symplectic basis of three-cycles A i: B\ with A4 fl = 8?. 
Here, the A4 cycles are the n three-spheres, and they project to cycles A4 surrounding 
the cut Ci = [x2i,X2i-i] in the x-plane. The Bi cycles are dual cycles which project in 
the x plane to the Bi paths The periods of Q are then given by 

ti = -i- 1 n, ^ = f a (3.80) 
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It is easy to see that these periods reduce to the periods ()3.63|) and ()3.75|) on the 
hyperelliptic curve ([3.74)1 . respectively. Let us consider again the case of the deformed 
conifold ()3.77|) . which is simpler since there is only one three-sphere. Let us compute 
the A-period over this three-sphere, which is an S 2 fibration over the cut [— ■JJI, \/P\i 
by first doing the integral over S 2 , and then doing the integral over the cut. Since 
v = y / p — x 2 — p 2 , where p 2 = y 2 + u 2 , the integral of Q over S 2 is simply 

1 f dydZ ^x~ 2 . (3.81) 



2tt Js 2 \/ ' p — x 2 — p 2 
Therefore, the A-period becomes 

1 



y(x)dx, (3.82) 



where y is now given by y 2 + x 2 = p. This is nothing but the A-period ()3.63|) (up to 
a redefinition y — > —iy). The general case is very similar, and one finally obtains that 
the special geometry ()3.80|) of the deformed Calabi-Yau geometry ()3.78|) is equivalent 
to the planar solution of the matrix model, given by the hyperelliptic curve 1)3.74)1 and 
the equations for the partial 't Hooft couplings (j3.64j) and the planar free energy (|3.75|) . 

The physical interpretation of this result is that there is an equivalence between an 
open topological string theory on the manifold X rcs , with N D-branes wrapping the 
n spheres obtained by blowup, and a closed topological string theory on the manifold 
X def , where the N D-branes have disappeared. Moreover, the 't Hooft couplings in 
the open string theory become geometric periods in the closed string theory. Since 
the open topological strings on X res are described by a matrix model, the fact that the 
planar solution reproduces very precisely the deformed geometry is important evidence 
for this interpretation. This duality relating an open and a closed string theory is an 
example of a geometric, or large N, transition. Notice that, as a consequence of this 
duality, the 't Hooft resummation of the matrix model corresponds to a closed string 
theory propagating on X^ e {. The master field controlling the planar limit (which is 
encoded in the planar resolvent, or equivalently in the quantity y(X)) leads to an 
algebraic equation that describes very precisely the target of the closed string theory 
dual. The large N transition between these two geometries was proposed in [15] . The 
fact that the open string side can be described by a matrix model was discovered in 



3.6 Extensions and applications 

The results derived above can be extended to more complicated Calabi-Yau back- 
grounds with branes PHIEHI- For example, one can consider ADE type geometries with 
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branes wrapping two-spheres [IHllIlj, and the string field theory description reduces to 
the ADE matrix models considered in |46j . In the one-matrix model described before, 
the master field is given by a hyperelliptic curve F(x, y) = which is then regarded as 
the Calabi-Yau manifold 

uv + F{x,y) = (3.83) 

in disguise. In some of the examples considered in [2BJ EH] , however, the master field is 
no longer described by a hyperelliptic curve, but involves a more complicated geometry. 
This geometry is the Calabi-Yau closed string background that is obtained by geometric 
transition from the open string background with branes. A detailed study of the more 
complicated master field geometries that arise in multimatrix models can be found in 

m- 

Another consequence of the result of Dijkgraaf and Vafa, together with the geometric 
transition of 15J, is that the Kodaira-Spencer theory of gravity j^j on the noncompact 
Calabi-Yau manifold ()3.78|) is equivalent to the 't Hooft resummation of the matrix 
model with potential W(x). For the simple example of the cubic potential, this was 
explicitly checked at genus one in jiS]. The formalism developed in [32] seems to be 
very appropriate to establish this equivalence in detail. 

As we mentioned in the introduction, the main application of the results of Dijkgraaf 
and Vafa has been the computation of effective superpotentials in supersymmetric 
gauge theories by using matrix model techniques. This is based on the fact 
that the resummation F (t) of the open string amplitudes is deeply related to the 
superpotential of the gauge theory which can be obtained from string backgrounds 
with branes. We refer the reader to [HI E2] for an exposition of these results. 

4 Type A topological strings, Chern-Simons theory 
and matrix models 

The conceptual structure of what we have seen in the B model is the following: first one 
shows, by using string field theory, that the target space description of open topological 
B strings reduces to a matrix model in certain backgrounds. Then one solves the model 
in the planar limit, and a geometry emerges which is interpreted as a closed string dual 
to the original open string theory. Both geometries are related by a large N transition. 
The first transition of this type was discovered in the context of topological A strings 
by Gopakumar and Vafa j3H|. What we will do here is to rederive their result by using 
the language and technology of matrix models. The key ingredient is the fact pointed 
out in that the partition function of Chern-Simons theory can be written in terms 
of a somewhat exotic matrix model. We will only focus on the matrix model aspects 
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of this correspondence. A detailed review of Chern-Simons theory and the geometric 
transition for the A model can be found in |54„. 



4.1 Solving the Chern-Simons matrix model 



S — — 

4:71 



The Chern-Simons action with gauge group G on a generic three-manifold M is defined 
by 

/ Tt(a AdA + ^-A A A A A) (4.1) 
Jm V 3 / 

Here, k is the coupling constant, and A is a G-gauge connection on the trivial bundle 

over M. We will consider Chern-Simons theory with gauge group G = U(N). As 

noticed in since the action (|4.1jl does not involve the metric, the resulting quantum 

theory is topological, at least formally. In particular, the partition function 

Z(M) = J [VA]e iS (4.2) 

should define a topological invariant of M. A detailed analysis shows that this is in fact 
the case, with an extra subtlety related to a choice of framing of the three-manifold. 

The partition function of Chern-Simons theory can be computed in a variety of ways. 
In [21] it was shown that in fact the theory is exactly solvable by using nonperturbative 
methods and the relation to the Wess-Zumino-Witten (WZW) model. In particular, 
the partition function of the U(N) theory on the three-sphere S 3 is given by 

z< - s3) = wtn^ £ ^M-ftV ■ "<4 (4 - 3) 

where the sum over w is a sum over the elements of the Weyl group W of U(N), e(w) is 
the signature of w, and p is the Weyl vector of SU(N). By using Weyl's denominator 
formula, 

<w)e w ^ u = Y[ 2 sinh (4.4) 

tueW a>0 

where a are positive roots, one finds 

where _ . 

2m 

»■ = *+»■ (4 - 6) 

It was found by Witten that open topological type A strings on T*S 3 (which is nothing 
but the deformed conifold geometry (|3.77jl ) in the presence of N D-branes wrapping 
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S 3 are in fact described by U(N) Chern-Simons theory on S 3 [74J. This is the type A 
model analog to the fact that open type B strings on the geometry described by (|3.4U|) 
are captured by a matrix model, and in both cases this is shown by using open string 
field theory. The free energy of Chern-Simons theory on S 3 has an expansion of the 
form (jl.2j) . with g s given in ()4.6|) . and the coefficients F 9t h, which can be computed 
by standard perturbation theory, have the interpretation of open string amplitudes on 
T*S 3 . 

The analogy between the A story and the B story can be taken even further, since it 
turns out that the partition function of Chern-Simons on S 3 , as well as on many other 
three-manifolds, can be represented as a matrix integral [SB]- In the case of S 3 most 
of the physical information in Z(S 3 ) can be obtained by other means, but for other 
three-manifolds like lens spaces and Seifert spaces, the matrix model representation is 
crucial in order to extract the coefficients F g ^ [SB]- The Chern-Simons matrix model 
on S 3 gives however a particularly clean way to derive the resummed free energies F g (t) 
and the geometry of the master field, and we will devote the rest of these lectures to 
presenting this analysis. 

In the case of S 3 the easiest way to derive the matrix model representation of the 
Chern-Simons partition function is through direct computation. Consider the following 
integral: 



N 

_ n^-E^n^sinh^)-. ,4.7) 

i=l i<j 

It can easily be seen that this reproduces the partition function of U(N) Chern-Simons 
theory on S 3 , given in (|4.5j) . and the derivation is left as an exercise. 

Exercise. Use the Weyl formula (|4.4jl to write (|4.7jl as a Gaussian integral, and 
show that it reproduces (|4.3|) . 

The measure factor in (|4.7|) 

n(2sinh^^) 2 (4.8) 

is not the standard Vandermonde determinant, although it reduces to it for small 
separations among the eigenvalues. In fact, for very small g s , the Gaussian potential in 
(14. 7|) will be very narrow, forcing the eigenvalues to be close to each other, and one can 
expand the sinh in (J4.8)) in power series. At leading order we find the usual Gaussian 
matrix model, while the corrections to it can be evaluated systematically by computing 
correlators in the Gaussian theory. In this way one obtains the perturbative expansion 
of Chern-Simons theory, see for details. 
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Here we will take a slightly different route in order to analyze the model. First of 
all, we want to write the above integral as a standard matrix integral with the usual 
Vandermonde discriminant. This can be achieved with the change of variables [HZ] 

exp(A + 1) = A i; (4.9) 

where t = Ng s , as usual. It is easy to see that the above integral becomes, up to a 
factor exp(— N 3 g s /2), 



1 r N N 

Zsw =M 11^ A '( A ) exp(-^(logA,) 2 /2^), (4.10) 
' J i=i i=i 

therefore we are considering the matrix model 

We will call this model the Stieltjes- Wigert matrix model, hence the subscript in ()4.10|) 
and ()4.11j) . This is because it can be exactly solved with the so-called Stieltjes- Wigert 
polynomials, as we will explain in a moment. 

Matrix integrals with logarithmic potentials are somewhat exotic, but have appeared 
before in connection with the Penner model [HI], with the c = 1 string at the self-dual 
radius jHHUj, and with the IP 1 model [Hi. We want to analyze now the saddle-point 
approximation to the matrix integral ()4.7j) . or equivalently to ()4.10j) . Since the model 
in ()4.10|) has the standard Vandermonde, we can use the techniques of section 2.2. 
Although the formulae there were obtained for a polynomial potential, some of them 
generalize to arbitrary polynomials. In particular, to obtain the resolvent ljq(p) we can 
use the formula ()2.60|) with 

W'(z) = (4.12) 

z 

Notice that this potential has a minimum at z = 1. We then expect a one-cut 
solution where the endpoints of the interval a(t), b(t) will satisfy a(0) = 6(0) = 1. In 
order to compute the integral ()2.60|) we deform the integration contour. In the case 
of polynomial potentials, we picked a residue at z = p and at infinity. Here, since the 
logarithm has a branch cut, we cannot push the contour to infinity. Instead, we deform 
the contour as indicated in Fig. |HJ we pick the pole at z = p, and then we surround 
the cut of the logarithm along the negative real axis and the singularity at z = with 
a small circle C e of radius e. This kind of situation is typical of the solution of matrix 
models with the character expansion The resulting integrals are: 

^ r dZ +/ dzl ° gZ ) (4.13) 



2*1 7-oo z(z -p)sj{z — a)(z — b) Jc € z(z - p)y / (z - a)(z - b) J 
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Figure 8: This shows the deformation of the contour needed to compute the planar 
resolvent of the Chern-Simons matrix integral. We pick a residue at z — p, and we 
have to encircle the singularity at the origin as well as the branch cut of the logarithm, 
which on the left hand side is represented by the dashed lines. 



Both are singular as e — > 0, but singularities cancel, and after some computations one 
finds for the resolvent: 



<*>o(p) 



1 

2tp 



log 



{y/ay/p -b - \/by/p - af 
(y/p — a — yjp — b) 2 p 2 



2tp\/ ab 



Aab 



2\/ab + a + b 
(4.1 

In order to satisfy the asymptotics (J2.57|) the second term must vanish, and the first 
one must go like 1/p. This implies 

2\/ab + a + b, 



Aab 



y/a~+Vb = 2e\ 



(4.15) 



and from here we obtain the positions of the endpoints of the cut a, b as a function of 
the 't Hooft parameter: 

3i 



a{t) = 2e 2t -e' + 2e^v / e F3 I, 
b(t) 



2e 2t - e*- 2e^\/e i - 1. 



(4.16) 



Notice that, for t = 0, a(0) = 6(0) = 1, as expected. The final expression for the 
resolvent is then: 



^o(p) = -— log 
tp 



1 + e"*p + ^(1 + e~ l p) 2 - Ap 
2p~ 



(4.17) 



and from here we can easily find the density of eigenvalues 



P(A) 



1 

TrtA 



tan 



'y/4X- (l + e-*A) 2 " 



1 + e-*A 



(4.18) 
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If we now define 

u(p) =t(l -pu (jp)) + ni 
we see that it solves the equation 



(4.19) 



v—u+t 



+ 1 = 



(4.20) 



where we put p = e t ~ v '. This was found in [2] by a similar analysis. The equation (|4.20j) 
is the analog of (|3.74jl in the case of polynomial matrix models, and can be regarded 
as an algebraic equation describing a noncompact Riemann surface. In fact, (J4.20j) 
is nothing but the mirror of the resolved conifold geometry (see for example |4TH ) , 
and t is the Kahler parameter of the geometry. This is of course in agreement with 
the result of 35J, who argued that the 't Hooft resummation of Chern-Simons theory 
leads to a closed string theory propagating on the resolved conifold. As in the B model 
that we analyzed before, the master field of the matrix model encodes the information 
about the target geometry of the closed string description, and provides evidence for 
the geometric transition relating T*S 3 and the resolved conifold geometry. 

As we mentioned before, the matrix model 1)4.11)1 can be solved exactly with a set 
of orthogonal polynomials called the Stieltjes-Wigert polynomials. The fact that the 
Chern-Simons matrix model is essentially equivalent to the Stieltjes-Wigert matrix 
model was pointed out by Tierz in |fi7j . The Stieltjes-Wigert polynomials are defined 
as follows 



2 i n X — Tli v{v — ri) o 1 



V 



(4.21) 



v=0 L J 



and satisfy the orthogonality condition ()2.87j) with 



i_( log 3 .)2 dx 
2tt 



(4.22) 



and 




where 




(4.23) 



' m [m] ! [n — m] ! 
The recursion coefficients appearing in ()2.92j) are in this case 



(4.24) 



r n = q 3n (q n -l), 



q h+n( q n+l +q n_ 1 y 
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The Stieltjes-Wigert ensemble can be regarded as a g-deformation (in the sense of 
quantum group theory) of the usual Gaussian ensemble. For example, as g s — > one 
has that [n] — > ng s , therefore 

K -> h%, (4.25) 

where is given in ()2.96|) . Also, one can easily check that the normalized vev of 
Tin M in this ensemble is given by 

(Tr R M) sw = e^g-dim, R, (4.26) 

where £(R) is the number of boxes of R, k r is a quantity defined by 

K R = £{R) + Y, X ^- 2i ) ( 4 - 27 ) 

i 

in terms of lengths of rows A; in R, and dim g R is the quantum dimension of the 
representation R 

dim,fi = n ' a ', (A ? ) ' I"' 28 * 
- LJ - \a ■ p\ 

where A is the highest weight associated to R. As g s — ■> 0, the vev (|4.26j) becomes 
just dimi?, the classical dimension of R, which is essentially the vev in the Gaussian 
ensemble (j2.8j) . 

Notice that, for this set of orthogonal polynomials, the expansion (j2.107|) is very 
simple since 

i? (O = e 44 «(l-e-*«), # 2s (£) = 0, s>0, 
s(0 = e*(l-2e*). (4.29) 

As we pointed out in section 2.3, i?o(0 an( i s (0 can De used to determine the endpoints 
of the cut in the resolvent through (|2.111j) . It is easy to see that ()4.29|) indeed lead to 
(I4.16|) . and that by using (|2.11(J|) one obtains (|4.18j) . In fact, it is well-known that the 
expression ()4.18|) is the density of zeroes of the Stieltjes-Wigert polynomials |4*H1 IT^j. 
We can now use the technology developed in section 2.3 to compute F g {t). Since 

Fes = F sw ~ ^t 3 + (4.30) 



the formula (j2.108j) gives 



F CS (t) = ^ " ^ " Li8(e-*) + C(3), (4.31) 



55 



where the polylogarithm of index j is defined by: 

oo n 
n=l 

The above result is in precise agreement with the result in [35] obtained by resumming 
the perturbative series. With some extra work we can also compute F^ s (t), for all 
g > 0, starting from (j2.112|) . We just have to compute /^(l) — /^(0), for p odd, 
where 

/(0 = (i - O0(e, *), <K£, *) = log + 4te. 



It is easy to see that 



and by using the expansion 



1 

7 + Em 



1-e- t ^ ' -"(*+!)! 
one gets 

0W(O,t) = fcl^tP. 

P 

Putting everything together, we find for g > 1 



B2gB 2g -2 B 2g , t . _ ^2g ^ 2 -2g 

2< 7 (2< ? -2)(2< ? -2)! 2g(2g-2)\ 3 " M ; 2g(2g - 2) 



Since the last piece is the free energy at genus g of the Gaussian model, we conclude 
that the Chern-Simons free energy at genus g is given by 

F > {t) - 2,3(23 -m 9 - 2)! + 2 j( 2 9 -2)! L ' 3 - 2 » (e > (4 ' 33) 

which agrees with the resummation of and also with the genus g closed string 
amplitude of type A topological strings on the resolved conifold (see [S3j for more 
details). 



4.2 Extensions 

We have seen that the matrix model reformulation of Chern-Simons theory provides 
an efficient way to obtain the master field geometry and to resum the perturbative 
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expansion. The result ()4.33|) can be derived as well from the perturbation series jHSElI, 
but the existence of a matrix model description of Chern-Simons theory turns out to be 
useful in other situations as well. For example, one can easily write a matrix integral 
for Chern-Simons theory for other gauge groups jSSj, arid the corresponding models 
have been analyzed in [37J. Moreover, the matrix representation of Chern-Simons 
partition functions can be extended to lens spaces and Seifert spaces, and provides a 
useful way to study perturbative expansions around nontrivial flat connections. The 
matrix models that describe these expansions have been studied in perturbation theory 
in [HTT^ E| and the saddle-point approximation to lens space matrix models has been 
studied in jSH]. There are as well multimatrix models describing A topological strings 
on some noncompact Calabi-Yau geometries j2] that can be studied by using saddle- 
point techniques [76], and it is possible as well to formulate the Chern-Simons partition 
function on S 3 in terms of a unitary model jHHj- However, all these matrix models are 
usually much harder to analyze than conventional ones, and more work is needed to 
understand their large N properties. 
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